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In this report we introduce the basic techniques (of the closed-time-path coarse-grained 
effective action) and ideas (scaling, coarse-graining and backreaction) behind the treat- 
ment of quantum processes in dynamical background spacetimes and fields. We show how 
they are useful for the construction of renormalization group (RG) theories for studying 
these nonequilibrium processes and discuss the underlying issues. Examples are drawn from 
quantum field processes in an inflationary universe, semiclassical cosmology and stochastic 
gravity. In Part I we begin by establishing a relation between scaling and inflation, and 
show how eternal inflation (where the scale factor of the universe grows exponentially) can 
be treated as static critical phenomena, while a 'slow-roll' or power-law inflation can be 
treated as dynamical critical phenomena. In Part II we introduce the key concepts in open 
systems and discuss the relation of coarse-graining and backreaction. We recount how the 
(in-out, or Schwinger-DeWitt) coarse-grained effective action (CGEA) devised by Hu and 
Zhang can be used to treat some aspects of the effects of the environment on the system. 
This is illustrated by the stochastic inflation model where quantum fluctuations appearing 
as noise backreact on the inflaton field. We show how RG techniques can be usefully applied 
to obtain the running of coupling constants in the inflaton field, followed by a discussion of 
the cosmological and theoretical implications. In Part III we present the Closed- Time-Path 
(CTP, in-in, or Schwinger-Keldysh) CGEA introduced by Hu and Sinha. We show how to 
calculate perturbatively the CTP CGEA for the model. We mention how it is useful 
for calculating the backreaction of environmental fields on the system field (e.g. light on 
heavy, fast on slow) or one sector of a field on another (e.g. high momentum modes on low, 
inhomogeneous modes on homogeneous), and problems in other areas of physics where this 
method can be usefully applied. This is followed by an introduction to the infiuence func- 
tional in the (Feynman- Vernon) formulation of quantum open systems, illustrated by the 
quantum Brownian motion models. We show its relation to the CTP CGEA, and indicate 
how to identify the noise and dissipation kernels therein. We derive the master and Langevin 
equations for interacting quantum fields, represented in the works of Lombardo and Mazz- 
itelli and indicate how they can be applied to the problem of coarse-graining, decoherence 
and structure formation in de Sitter universe. We perform a nonperturbative evaluation 
of the CTP CGEA and show how to derive the renormalization group equations under an 
adiabatic approximation adopted for the modes by Dalvit and Mazzitelli. We assert that 
this approximation is incomplete as the effect of noise is suppressed. We then discuss why 
noise is expected in the RG equations for nonequilibrium processes. In Part IV, following 
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Lombardo and Mazzitelli, we use the RG equations to derive the Einstein-Langevin equation 
in stochastic semiclassical gravity. As an example, we calculate the quantum correction to 
the Newtonian potential. We end with a discussion on why a stochastic component of RG 
equations is expected for nonequilibrium processes. 



I. AIM AND SCOPE 



We discuss how the concepts of open systems and the techniques of the coarse-grained effective action 
and influence functionals can be applied to nonequilibrium quantum processes in the early universe and in 
semiclassical gravity, leading to renormalization group (RG) theories for the description of the interaction 
dynamics of these theories. 

Our wish is that by examining a sample class of problems of fully dynamic nature -as different from 
equilibrium (finite temperature) or near-equilibrium (linear response theory) - we can lay out the issues and 
approaches useful for the construction of a RG theory for nonequilibrium (NEq) processes involving quantum 
fields. 

Description of phase transitions involves a scale that measures the behavior of the order parameter field 
in the critical region. The energetics of the system is characterized by its quantum dynamical and statistical 
mechanical properties. At the heart of NEq statistical mechanics is the interplay of the dynamical scales of 
the system (from the time-dependent order-parameter field) and some background (e.g., all physical processes 
in an expanding universe are measured against the time-dependent metric function such as the scale factor 
a{t)). In equilibrium treatment, this is usually captured by a finite temperature effective potential. But 
since the order parameter field is generally time dependent, and there may not be a thermal equilibrium 
environment present in these dynamical processes, one should really be working with an effective action or a 
free energy density functional. Unlike scattering problems commonly found in particle physics where one can 
determine the transition amplitude between the in and out states based on the in-out or Schwinger-DeWitt 
effective actions, in evolutionary problems frequently encountered in statistical mechanics, the development 
of the expectation value (of an operator associated with some physical variable, such as the energy momentum 
tensor) need be obtained from the in-in or Schwinger-Keldysh effective actions. When open system concepts 
like coarse-graining and backreaction of the system and environment are applied, the CTP CGEA or the 
influence action are more appropriate. (For the development and application of these ideas applied to 
problems in gravitation and cosmology, see, e.g., ||l|-§l). 

This paper is in the nature of a report rather than a review - in that we will present or develop only 
those works which are useful for the construction of the conceptual and technical frameworks to treat this 
broad class of problems, taking specific examples from semiclassical gravity and inflationary cosmology as 
illustrations. 

RG concepts and techniques have been used in other areas of gravitation and cosmology. We mention 
some representative works; 

1) RG in gravitational collapse: The use of universality and scaling ideas in classical gravitational collapse 
first discussed by Choptuik ||^ has grown since then into an interesting area of classical graviational research. 
Fractal structure and scaling laws in a self-gravitating gas have been investigated by de Vega, Sanchez and 
Combes |lO) 

2) RG in quantum field theory in curved spacetime: notable work since the 80's by Calzetta, Hu, O'Connor, 
Jack, Parker, Toms and others as well as the Tomsk group can be found in the book of Buchbinder, Odintsov 
and Shapiro [0. For more recent works see 

El and 

3) Scaling in quantum gravity has been pursued by Ambjorn in a simplicial gravity approach and by 
Antoniadis, Mazur and Mottola |15|, to mention just a few notable avenues of inquiry. 

The organization of this paper is as follows: In Sec. 2 we first give a descriptive summary of the nature of 
the problems encountered in phase transitions in the early universe, focussing on the inflationary cosmology. 
We distinguish the case of eternal inflation from that of slow roll and indicate why they can be viewed 
as static and dynamic critical phenomena respectively. In Sec. 3 we begin a discussion of the relation of 
scaling, coarse-graining and backreaction with the example of stochastic inflation, thus bringing out the basic 
concepts of open systems. In Sec. 4 we introduce the ('in-out') coarse-grained effective action (CGEA) for 
a A(/)^ field to incorporate the backreaction effect of the short wavelength modes viewed as the environment, 
on the long wavelength modes of the system. In Sec. 5 we perform a rescaling of the modes and the field 
in the spirit of RG transformations and derive the correspondng RG equations for the coupling constants 
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of the system field. We briefly discuss the theoretical and cosmological implications. In Sec. 6, we present 
the 'in-in' or closed-timc-path (CTP) CGEA and show how it is useful for the derivation of real and causal 
dynamical equations for expectation values of operators of the system fields with backreaction from the 
environment field. We carry out a perturbative evaluation of the CTP CGEA and show how it could be 
useful for the consideration of backreaction of environment field (modes) on the system field (modes). In 
Sec. 7 we introduce the influence functional formulation via the quantum Brownian motion model and show 
its relation with the CTP CGEA. In this process we obtain the master and semiclassical Langevin equations 
for interacting quantum fields. We show how the noise kernels can be identified in this equation and how 
their behavior can be used as a measure of decoherence. In Sec. 8 we carry out a nonpcrturbative evaluation 
of the CTP CGEA, and derive the RG equations under different approximations. We indicate how our real- 
time CTP Coarse-Grained Effective Action (CGEA) approach differs from the Euclidean Averaged Effective 
Action (EAEA) approach. In Sec. 9 we use RG theory to derive the Einstein-Langevin equation in stochastic 
semiclassical gravity. In Sec. 10 we show how the RG equations change the Newtonian potential. In Sec. 
11, following a short summary, we discuss the salient features of RG theory for nonequilibrium processes and 
argue why there is a stochastic component in the RG equations for such systems. 

Those readers who only wish to learn the methodology which underlies nonequilibrium quantum field 
processes, but have no special interest in the specifics of such processes in gravitation and cosmology, can 
do just read Sections 4, 6, 8 and 11, from which they will be able to apply the methods to their own areas 
of research (e.g., nuclear/particle, atomic/optical physics). 



PART ONE: Scaling and Inflation 

II. PHASE TRANSITIONS IN THE EARLY UNIVERSE 

We use phase transitions in the early universe (from the Planck time to the GUT time, or even at 
the later electroweak era) as working examples to illustrate the basic issues and methods involved. This 
is because phase transitions in the early universe are mediated by and involve many NEq processes of 
fundamental interest (e.g., nucleation, spinodal decomposition, particle creation, decoherence) that lead 
to many important physical consequences such as entropy generation and structure formation. We will 
only dwell on those aspects which illustrate the applications of the RG theory and suggest its extension 
to dynamical processes. The development of a RG theory for NEq processes here - from the construction 
of the CGEA to the derivation of the RG equations — is based on the concepts (scaling, coarse-graining, 
backreaction) and techniques (CTPEA and influence functional) of nonequilibrium mechanics and quantum 
field theory. These ideas and techniques are generally applicable also to problems outside of gravitation 
and cosmology, such as the quark-giuon plasma and atoms interacting with a Bose-Einstein condensate. On 
the other hand this parallel presentation of RG theory in the context of NEq processes may be helpful to 
cosmologists who wish to find a more solid theoretical anchor for the discussion of phase transitions in the 
early universe. 

A. Effective Action for Dynamic Order-Parameter Fields 

As noted earlier knowledge of the exact form of the CTP CGEA holds the key to a complete description 
of a phase transition. One can deduce not only the qualitative features (first or second order) but also 
the quantitative details (mechanisms and processes) . Therefore the construction of the effective action for 
the order parameter field in different cosmological spacetimes is usually the necessary first step towards a 
description of phase transitions in the very early universe. 

For the purpose of illustration it is useful to establish the connection with ideas of RG theory at the 
outset. The two themes established in QJ^, which we follow here are: 1) Eternal infiation can be described 
equivalently as an exponential scale transformation, thus rendering this special class of dynamics as effectively 
static. 2) The class of 'slow-roll' infiation can be treated as a dynamical perturbation off the effectively static 
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class of exponential inflation and be understood as a dynamical critical phenomenon in cosmology, q 

Let us concentrate on situations where the order parameter field changes either with space or time. (A 
familiar example in condensed matter physics is anisotropic superconductivity where one can use a gradient 
expansion in the Landau-Ginzburg- Wilson effective potential to account for the differences coming from 
the next-to-nearest neighbor interactions.) For cosmological problems, it is the time-dependence of the 
background field which one needs to deal with. |^ Thus for a realistic description of many inflationary 
transitions one needs to treat the case of a dynamical field and a nonflat or even quasi-static potential. 
The form of the potential and the metric of the background spacetime together determine the behavior of 
the scalar field in the Laplace-Beltrami equation, but the field in turn provides the source of the Einstein 
equation which determines the behavior of the background spacetime metric. Hence they ought to be solved 
self-consistently. (One usually considers only the homogeneous mode of the scalar field for the dynamics of 
inflation and the inhomogeneous modes of quantum fluctuations for processes like structure formation.) 

At the classical level, the wave equation for the background scalar field (assumed homogeneous) with 
self-interaction potential V{(j)) in a spatially-flat Robertson- Walker (RW) spacetime is given by 

4> + 3H{t)4> + V'{<t)) ^0 (III) 

and the Einstein equations read 

H + 3H^ = 87rGV{(j)), H = -AnG^^ (II.2) 

where H{t) = a/a is the Hubble rate, a dot denotes derivative with respect to cosmic time t, and a prime 
denotes a derivative taken with respect to its argument. A trivial but important solution to these equations 
is obtained by assuming that V((j)) = Vq — constant, (j) = = constant and H = Hq = constant, which 
is the de Sitter universe a = e^°* with a constant field. A less trivial but useful solution is the so-called 
'power-law' inflation models [2^ , with an exponential potential and a slowly- varying inflaton field. One 
can carry out a derivative expansion of the background field to obtain a quasi-local effective action for the 
description of such classes of spacetimes and fields. (An example of 'slow-roll' transition.) 



'^A comment on the meaning of the term 'critical dynamics' as used in the context of cosmological phase transitions 
is in order. By it we refer to studies of phase transitions mediated by a time-dependent order parameter field in 
contrast to static critical phenomena where the order-parameter field is constant in time. We are using this term in 
a general sense, not necessarily referring to the specific conditions of critical phenomena as discussed in condensed 
matter systems [^. For example, critical phenomena usually deals with the change of the order parameter field 
near the critical point as a function of temperature. In cosmology, temperature T is a parameter usually (e.g. under 
the assumption of adiabatic expansion) tied in with the scale factor a and does not play the same role as in critical 
phenomena. In the new inflationary scenario, the critical temperature Tc is defined as the temperature at which a 
global ground state (the true vacuum) first appears. The stage when vacuum energy begins to dominate and inflation 
starts is the beginning of the phase transition. The stage when the system begins to enter the true vacuum and reheat 
can be regarded for practical purposes as the end of the phase transition. Throughout the process of inflation the 
system is in a 'critical' state. The progression of a cosmological phase transition is measured not by temperature, but 
by the change of fleld configurations in time. Criticality corresponds to the physical condition that the correlation 
length ^ = rn~^j: ^ oo, or m^^j = ci^'F/d(^^ 1^=0 (which may or may not be possible). In critical dynamics studies 
of condensed matter systems one usually analyses the time-dependent Landau-Ginzburg equation, with a noise term 
representing the effect of a thermal bath, and studies how the system (order-parameter field) settles into equilibrium 
as it approaches the critical point. We are not concerned with the corresponding cosmological problem here. An 
attempt to describe this aspect of the inflationary transition was made in See also p^ . 

^Strictly speaking, phase transition studies usually carried out assuming a constant field in the de Sitter universe are 
unrealistic, in that they only address the situation after the universe has entered the inflationary stage and inflates 
indeflnitely. This model cannot be used to answer questions raised concerning the likelihood that the universe will 
still inflate if it had started from a more general, less symmetric initial state, such as the mixmaster universe |]l9|-pl[ . 
Nor can one use this model to study the actual process of phase transition (e.g., slow roll-over), and the problem 
of exit (graceful or not) to the 'true' Friedmann phase. To do this, as is well-known, one usually assumes that the 
potential is not exactly flat, but has a downward slope which enables the inflaton fleld to gradually (so as to give 
sufficient inflation) settle into a global ground state. The cosmological solution is, of course, no longer a de Sitter 
universe. 
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B. Inflation as Scaling: Static Critical Phenomena 



This idea arose from the work of Hu and Zhang Q on coarse-graining and backreaction in stochastic 
inflation. In trying to compare the inflationary universe with phase transitions in the Landau-Ginzburg 
model, using a field as example, they realized that the exponential expansion of the scale factor can be 
viewed as the system undergoing a Kadanoff-Migdal scale transformation ||2^ (this is explained in detail in 
Sec. 3.B). In other words, time in this case plays the role of a scaling parameter. It does not have to be 
viewed as a dynamical parameter. Thus for this special class of expansion, the dynamics of spacetime can 
be replaced equivalently by a scaling transformation. In so doing one renders eternal inflation into a static 
setting. By contrast, the larger class of power-law expansion a = f does not possess this scaling property. A 
useful parameter which marks the difference between these two classes of dynamics is C = \H\/H'^ = a/d^, 
where a = Ina, which can be regarded as a 'nonadiabaticity parameter' of dynamics: the de Sitter exponential 
behavior with C = is 'static', the slow-roll with small C << 1 is 'adiabatic', while the RW low-power-law 
with C ~ 1 is 'nonadiabatic'. 



C. Quasi-Static Field, 'Slow-Roll' as Dynamical Critical Phenomena 

The effective potential V{(l)) gives a well defined description of phase transition only for a constant back- 
ground (order-parameter) field. If the order-parameter field is dynamic, the effective potential is ill- defined 
and a host of problems will arise. Indeed, the very meaning of phase transition can become questionable. 
This is because as the field changes the effective action functional changes, and the location of the min- 
ima changes also. The notion of symmetry breaking and restoration is meaningful only when there exist 
well-defined global and local minima which do not change much in the time scale of the phase transition. 
Changing background fields will also engender particle creation, which affects the nature and energetics of 
phase transition as well. Therefore, in the context of phase transitions involving dynamic fields, short of 
creating a new framework, one can at best discuss the problem in a perturbative sense, where the back- 
ground field is nearly constant (quasi-static), so that an effective quasi-potential can still be defined |p^ , p5| . 
An effective Lagrangian for a slowly-varying background field can be obtained by carrying out a quasilocal 
expansion in derivatives of the field, the leading term being the effective potential p6t . 

C = C{cb,d^cf>,d^d,cb,...) (II.3) 

One can use this method to derive effective quasi-potentials for scalar fields in flat space (for an example 
of its application to electroweak finite temperature transition see, e.g., Moss et al |2^) or (in conformal 
time) for the conformally-fiat Robertson- Walker spacetimes. This is useful for studying cosmological phase 
transitions where the background spacetime changes only gradually, as in the Friedmann (low-power law) 
solutions a = n,7 < 1. (For a description see [^.) However, for the inflationary universe where the scale 
factor undergoes rapid expansion following either an exponential a — e^* or a high power-law behavior, the 
quasilocal expansion which assumes that the background field varies slowly is usually inadequate. 

Using the conceptual framework introduced above, one can understand why the particular subclass of 
high-power-law expansion associated with an exponential potential can hence be viewed as quasi-static. It 
is in this context that one can introduce the quasi- 'static' approximation to derive the effective action for 
scalar fields to depict this more realistic 'slow-roll' inflation, now carried out as a quasilocal perturbation 
from the de Sitter space, which is viewed as effectively static. For slowly-varying background fields one can 
use the method of derivative expansion to derive the quasilocal effective Lagrangian. Usually this makes 
sense only for static (or conformally-static, like the RW) spacetimes. However, one can view the special class 
of exponential expansion as effectively 'static'. This can be understood with the ideas of 'dynamical finite 
size effect' and implemented by treating inflation as 'scaling' transformations [Q. The 'slow roll-over' 
type of phase transition used in many inflationary models can be viewed as a quasi-static case, and derived 
as a dynamic perturbation from the de Sitter universe This view reveals a close analogy of this case 
with dynamical critical phenomena where the scaling parameter s plays the role of inflation and the time 
parameter t measures how much the system departs from the exponential expansion solution. The main 
points proposed in |^ can be summarized by the following schematic diagram: 

Constant Field in Static or SCALING Exponential Expansion a = e^* 
Conformally Static Spacetimes > 'Eternal Inflation' 
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(Finite Size Effect) 



(Dynamical Finite Size Effect) 



Quasilocal Approximation 



Derivative Expansion 



High Power-Law expansion a = f 
^> 'Slow Roll-Over' 



Slowly Varying Field in 
RW Universe (Low Power-Law) 

The inadequacy of the finite temperature effective potential for the description of phase transitions in the 
early universe was what motivated us to look for more general methods useful for dynamic and nonequilibrium 
processes, especially involving quantum fields. We will trace out this pathway with two themes, one for the 
exposition of methods to treat NEq processes, and the other for the development of RG theories for such 
processes. This will lead us from the effective action to the infiuence functional methods, and for RG theory, 
from CTP CGEA to RG equations. 

It is better to start with a physical example to motivate, so we first discuss a simple conceptual point 
which will enable us to see inflation in the light of scaling. |^ 

III. COARSE-GRAINING, SCALING AND INFLATION 
A. Inflation 

Consider a massive (m), self-coupled (A) scalar field (f>{x,t) in a background spacetime with action 



SI 



4! 



{111.1) 



where □ 



minimal and conformal coupling of the field with the background spacetime with scalar curvature R. We 
will consider Robertson- Walker (RW) and de Sitter (DS) spacetimes with line element 



g-J^ is the Laplace-Beltrami operator and ^ = 0, ^ correspond respectively to 



ds^ = dt^ - a^dn^ 



(IIL2) 



where dil^ is that of the 3-space {dft^ = dx^ for spatially fiat). We restrict our attention here only to the 
spatially homogeneous and isotropic RW, or the spacetime- homogeneous and isotropic DS cases, where a 
single isotropic scaling parameter is applicable. For the RW universe 

a{t) = F, 

where 7 = 1/2 for a radiation-dominated source and 7 = 2/3 for a pressureless dust source. For a de Sitter 
universe (in the spatially-fiat RW coordinatization) 



a{t) — exp(Ht). 



(IIL3) 



Spatially homogeneous but anisotropic spacetimes are more complicated as they require different scaling 
parameters in different directions. 

In a spatially-homogeneous spacetime, the scalar wave function ^{x, t) 



Sec. 3, 4, 5 are excerpted from Lectures I and II of 
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(IIL4) 



is the product of a function of time (j)k{t) and a function of space Uk{x)[= e^^'^ for spatially- flat, and 
= YimniX: ^ t 4') i ^ = {l,m, n) for spatially-closed cases.] The wave equation for the time-dependent amplitude 
function (pkit) associated with the kth mode of a scalar field with self-interaction potential V{^) in 

a RW spatially-flat spacetime is given by 

^ + mt)^ + + + mcbk + vM = (III.5) 

where H{t) = d/a — a , or a — In a = J Hdt. (Here a dot denotes derivative with respect to cosmic time t). 



The solution to III. 5 depends on V {(j)) — dV/dip. For inflation we are interested in cases where i) 
V'{4>) << 3H(j) . This would correspond to considering only the relatively flat portion of the potential, 
where inflation takes place for an extended period of time (usually assumed Aa — J^^ Hdt > 68 e-folding 
time to engender sufficient entropy corresponding to the observed universe) . In new inflation, this corresponds 
to assuming a flat plateau; and in chaotic inflation, a gradual downward slope. In the realistic de Sitter case, 
the order parameter field (j) changes very slowly in this inflation regime (slow roll-over), its rate controlled 
by 3Hdip/dt, which can be regarded as a viscous term in the dynamics. One can safely also assume that ii) 
(p << 3H(f> in this regime. (By contrast the (f) term dominates the dynamics in the reheating regime). A 
further simplification is to assume that (j) = constant in this regime throughout, corresponding physically to 
eternal inflation, i.e., a = e^* for all time. This idealized case is of course unphysical in that the universe will 
never roll-over (in new and chaotic inflation), or tunnel (in old inflation) to the true vacuum. Nevertheless 
it captures the salient features of inflation. We shall consider this idealized case in detail here in the light of 
static scaling. In summary one can distinguish the following cases: 

a) Minkowski field theory (static) H = 

b) Eternal inflation (stationary) H = constant, = 

c) Slow-roll inflation (quasi-stationary) H ~ constant,!/) ^ 0, = 

d) Realistic inflation (dynamic) H ~ constant, (j) 0, cj) 



In both cases a) and b), the term SHcf) in III. 5 is ignored, while in cases b) and c) the term (f> is ignored. 

We now discuss the effect of (eternal) inflation on the scalar field. We will show that a scalar field with 
At/)'* self-interaction in an isotropically expanding spacetime undergoing eternal inflation beliaves exactly like 
the same field in a static spacetime undergoing a scaling transformation. The physics of inflation (we refer 
specifically to the initial stage only, but not the subsequent slow roll-over nor the reheating stages) can thus 
be understood in terms of scaling completely. 



B. Scaling 



To fix ideas, consider the spatially-flat RW metric Eq. III.2 with a = constant. This is just the Minkowski 
spacetime. Let us consider an ordered sequence of such static hyperspaces (foliation) with scales ao,ai,a2, 
etc, parameterized by tn — ta+nAt, rt = 0, 1, 2, . . . . These spacetimes have the same geometry and topology 
but differ only in the physical scale in space. One can always redefine the physical scale length x'. 



(n) 



to render them equivalent. If each copy has scale length magnified by a fixed factor e^^* over the previous 
one in the sequence, i.e. a„+i/a„ = e^'^*, we get exactly the physical picture as in an eternal inflation. 
After n-iterations i.e. a„/ao — e"*^^'^*\ or, in terms of a continuous parameter a{t) = ape^*. It is important 
to recognize here that t can be any real parameter not necessarily describing the dynamics. Time in this 
depiction plays the role of a scaling parameter, and dynamics is nothing other than scaling. 

To see this simple point in another light, let us describe a totally different physical situation (which we 
will see at the end is exactly identical) where scaling plays the dominant role but one would otherwise not 
be inclined to associate with the physics of inflationary cosmology. This is the Kadanoff-Migdal transfor- 
mation ||2^,|2^ widely applied to condensed matter systems for the study of critical behavior. Consider a 
3-dimensional cubic lattice with lattice length L, and an order parameter field (j) describing the magnetization 
of the system. The Ginsburg-Landau- Wilson Hamiltonian for a Heisenberg magnet bears similarity with the 
A(/)* theory in a well-known manner | |3C[ |. 

The Kadanoff-Migdal transform is an artificial procedure for relating the microscopic and macroscopic 
properties of a system based on the existence of scaling properties in the system in the infrared limit. It 
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involves taking a certain number of fixed spin sites (e.g. 4 in 2 dimensions, 8 in 3 dimensions) and replacing 
them by one block spin with adjusted (e.g. doubled bond strength) interaction strength between nearest 
neighbors. Carrying out this transformation n times leads to a coarse-grained system. If the system (ferro- 
magnet) manifests scaling properties near the critical point, as had been observed before the renormalization 
group theory was invented for its description, then the resulting rescaled macroscopic system should preserve 
the same properties as the original microscopic system. The appearance of long range order near the critical 
point makes this procedure a viable one, permitting enormous simplification while preserving the salient 
features of the system. 

Let us examine a two-dimensional example of this process, sometimes also called decimation. Denote by a 
the original length and a„ the scale length after carrying out the K — M transformation n times. If 4-spins 
are combined into 1 block spin in 2 dimensions, oi = 2ao, and a„/ao = 2". This gives the same result as 
inflation, except that the scaling factor s = a„+i/a„ here is 2, instead of s = e^^^^K Defining 

, ^ ^ ^ e^ (III.6) 



cr for inflation is HAt , a constant, while for decimation in the above example is In 2. We shall see that it 
is precisely this scaling property in inflation which imparts all the distinct and desirable physical features, 
from the scale-invariant Zel'dovich-Harrison perturbation spectrum to the Hawking effect. 

While the above real-space scaling explains well the physical idea, the dual transformation in momentum 
space is easier to calculate. This is the original Wilson-Fisher RG transformation |^ In momentum 

space, the lattice spacing acts like an ultraviolet cut off A of wave numbers A: (A = oo for continuous systems). 
The block spin transformation corresponds to a) eliminating the higher wave number modes (e.g. k > A/s) 
and b) rescaling (e.g., k —> k' = ks). The real parameter s > acts in a) like a coarse-graining (cut 
parameter) and in b) like a rescaling parameter. In our example above, s = 2. In every iteration, wavemodes 
with k > A/2 are integrated out first, then the remaining modes from < fc < A/2 are rescaled such that 
the new fc' space covers the full range with the same UV cutoff. For definiteness, we will refer to these two 
steps as coarse-graining (separation) and rescaling. Together these two steps constitute the renormalization 
group transformation. 

Thus in this example, at any iteration, wave numbers A/s < fc < A are being integrated away, and only 
the lower wave mode sector is kept. For the long wavelength, low fc cutoff, we assume the horizon size H~^, 
because those with lower fc having wavelengths greater than the horizon size are of no physical significance 
to the system, at least during the inflationary stage. (Some of these wavelength will reenter the horizon in 
the reheating phase and influence later events.) In the language of subdynamics one can call the former the 
environment (irrelevant sector), and the latter the system (relevant sector). The division of these two sectors 
changes with the order of iteration (in the language of Kadanoff-Migdal transformation) or time-evolution 
(for inflationary cosmology) of the system. The coupling of the environment with the system is determined 
by V{(l)). In addition, for a self-interacting field there are low-low and also high-high mode couplings. A 
useful way to keep track of these interactions, in particular, is by way of the so-called coarse-grained effective 
action which we will discuss in the next section. 



C. Coarse-graining and Stochastic Inflation 



In the above discussion we have shown how one can transcribe the dynamical process of inflation to the 
static transformation of scaling. We also showed how the evolution of the universe during inflation can be 
viewed as a succession of coarse-grainings, with the higher wave number modes cast away and their effect on 
the remaining (low fc )modes accounted for by a coarse-grained effective action. The resulting RG equations 
describe the flow of the coupling constants. Using this equivalence, we can conceptually (and technically) 
replace the inflationary expansion process in time by a running of the field towards the infrared regime. 

Let us compare the coarse-graining scheme used in the critical phenomena example above (GP, or Case A) 
with that of an important class of inflationary cosmologies called stochastic inflation (SI, or Case B), because 
there are close similarities and basic differences. Starobinsky's scheme is to divide a free scalar field into 
two parts: the system contains long- wavelength modes outside the horizon (p = k/a < eH,£ = being 
the physical wavelength and e < 1 a small parameter), which are treated classically; and the bath contains 
quantum fluctuations of wavelengths £ shorter than the horizon. Thus the bath in stochastic inflation is 
obtained from the full spectrum by a window function 6{k — eaH), and changes at each moment. Many 
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authors have claimed that these high k modes with a dynamic cutoff which obey the wave equation (Iff. 5) 
(with first order or with both first and second order time derivative terms) give rise to a white noise. They 
then proceed to set up and solve a Langevin equation for the classical field (f) driven by the de Sitter dynamics 
and with this white noise source. Thus the name stochastic inflation. Q 

The physical setup of Starobinsky's stochastic inflation differs from the example we used above to illus- 
trate the relation of scaling and inflation in how the system is defined. Technically there exists a simple 
transformation between coarse-graining in critical phenomena (CP, Case A) and that in stochastic inflation 
(SI, Case B). Let us explain: 

Compare the composition of the system (S) and the environment (E) at the beginning and end of inflation 
for critical phenomena (Case A) and for stochastic inflation (Case B). For Case B, at the start to, assume 
a{to) = 1 and e = 1 for simplicity. At any later time the system consists of all modes k < Ha or i > H~^. 
As the universe inflates a = e^*, more k modes are being red-shifted out of the horizon into the system. 
Thus S of stochastic inflation increases (from the k region HqK at to to HK at to GK at , see 
Fig. 1). The reason why in the stochastic inflationary program the system is taken to consist of modes with 
wavelengths greater than the horizon H^^ in the inflation regime is because these are the modes which would 
reenter the RW horizon in the radiation-dominated RW phase; and only those higher {k > kc) modes (here 
treated as environment ) which reenter the horizon before galaxy formation can physically influence fee in 
the RW phase. So the system in stochastic inflation is really designed for and chosen by the physics in the 
later radiation-dominated RW regime, rather than during the inflationary regime. Note that modes with 
wavelengths greater than the horizon size at the onset of inflation have not yet reentered today's horizon, 
because in, say, the new inflationary scenario, today's observable universe is believed to be only a small 
portion (Ih) of a much larger bubble inflated from a fluctuation domain (Id) at tQ. The mode kn exits 
the horizon at time t^ . The wave mode kn of physical wavelength ^//(today's Hubble radius) which is 
just entering the RW horizon today [and which disappeared outside the de Sitter horizon earlier than any 
observable physical scale (HRhHyiow in Fig. 1)] forms the upper limit of ^pjiyg- Of course more wave 
modes with £ > £h will be entering into the RW horizon tomorrow. These modes have not yet come into 
contact with our universe so they should be considered as unphysical by today's observers. However, they 
had interacted previously in the de Sitter phase at times t < t^ with all other k modes < kng within the de 
Sitter horizon. In the inflationary phase, it is this interaction which one should address. 

By contrast, in critical phenomema (Case A) depicted earlier, both the system and the bath are contained 
in the horizon throughout the period of inflation from tg to ti (the OLHi portion in Fig. 1). is the 

lower (IR) cutoff. The partition is defined by the incremental fraction s ~ e^^^ = a{t + At)/a{t) of k 
modes slipping into the bath at each moment Ai(or each iteration) by coarse-graining. For example, one 
possible partition fc^ (point L in Fig.l) is determined by the last wave mode which left the horizon at the 
end of inflation at ti. All modes with k > ki, never made it out and always stayed inside of the DS horizon. 
These will be the common "bath" to all systems at all times in the DS phase. Another choice more geared 
towards the galaxy formation problem is defined by the k modes of interest (e.g. galaxy scale kc) and made 
up of A; < kc. The environment then contains modes k > ko, including the common portion k > kL- In 
this designati on we make sure that: 1) The system consists of all scales within today's Hubble radius (RW 
horizon). 2) Their interactions in the inflationary phase with the higher and lower modes are incorporated. 
3) The environment is always within the de Sitter horizon. We are interested in two effects: 1) The effect 
of the higher k modes field 0> on the system field made up of lower k modes 0<, and 2) How the field 
parameters (such as mass m and coupling constants A) change in this process of inflation (or scaling) . These 
are answered respectively by 1) deriving the effective action and the associated effective equations for (j)^ , 
and 2) deriving the renormalization group equation for these field parameters. 

Despite the difference in the aims of the two programs the coarse-graining technique used in both cases 
are related in a simple manner. By replacing s — > s~^, or equivalently (At) — > —(At) in Case A, one gets the 
procedure in Case B for stochastic inflation. This amounts to changing inflation to deflation (time- reversal), 
or interchanging system and bath. One can also view the difference as representing a passive versus an active 



^Authors of Pj6|j32) question how a white noise can be deduced by this dynamic partitioning of a free field. They 
suggested instead how colored noise can be generated from coarse-graining an interacting field without assuming any 
dynamic partitioning (see Section VII B). This provides a more general scheme for stochastic inflation, where one 
can discuss issues like non-Gaussian galaxy distribution from vacuum fluctuations. 
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view. By this we mean inflating the background (spacetime) has, for a local observer in a scale-invariant 
system, the same effect as coarse-graining the system in a fixed background. 

A figurative description which encapsulates the main ideas in this comparison is the effect of a zoom 
lens. The physics of a scale-independent system under inflation is equivalent to holding the system static 
and viewing it through a zoom lens. Here zooming (scaling) replaces dynamics (inflation). Features of an 
evolving system undergoing inflation (the passive interpretation) can be depicted equally by the magnifying 
view through a zoom lens (scaling) of a fixed system (the active interpretation). Obviously this relation 
holds only if the structure is homogeneous and the lens has uniform magnification. (This corresponds to 
homogeneity in space and uniform scaling.) 

In the next section we present an effective action formalism for carrying out the coarse-graining. We shall 
keep the formalism general and non-specific to the system-environment stipulation so that it can be applied 
to different situations in a variety of problems. 



PART TWO: Coarse- Graining and Renormalization Group 

IV. COARSE-GRAINED EFFECTIVE ACTION 

A. Coarse Graining and Backreaction 

In quantum field theory, when the background field or spacetime dynamics follows some simple discernable 
(e.g., classical) behavior one commonly performs what is called a background field decomposition, 

on a (e.g., self-interacting) scalar field 4>. One then calculates the effect of the fluctuation field 0^ on the 
background field 4>c (as well as the background geometry g^i, when necessary) using a loop expansion in the 
in-out effective action. For many problems such as those encountered in statistical mechanics where one is 
interested in the detailed behavior of only a part of the overall system (call it the system) interacting with its 
surrounding (call it the environment), one can decompose the field describing the overall system <f> according 
to 

= i's + 4>E 

where (j)s denotes the system field and 4>e the environment field. One can calculate the backreaction of the 
latter on the former with the 'coarse-grained effective action' (CGEA) first introduced by Hu and Zhang Q], 
where one does not integrate over the quantum fluctuations, as in the conventional definition (Schwinger- 
DeWitt) of effective action in orders of ?i, but rather with respect to a small parameter e which can generally 
be expressed as the ratio of two time, length, mass or momentum scales characteristic of the system and 
environment discrepancy (examples are slow-fast time scale, long-short length scale, light-heavy mass scale 
or soft-hard momentum scale). In the stochastic inflation scheme one regards the system field as 
containing only the lower modes and the environmental field as containing the higher modes with the division 
provided by the event horizon in de Sitter spacetime. In post- inflationary reheating, the quantum (in- 
homogeneous) fluctuation fields are parametrically amplified by the time-dependent (homogeneous) inflaton 
(background) field, resulting in particle creation. Its backreaction effect is to damp away the oscillations of 
the inflaton field and reheat the universe. (See, e.g., |3^,Q). Similar processes of particle creation at the 
Planck time leading to damping of anisotropy of the background spacetime were earlier treated by means 
of a closed-time-path (CTP), or in- in effective action. (See, e.g., In quantum cosmology one usually 
makes the so-called 'minisuperspace approximation' of quantizing only the homogeneous universes and ig- 
noring the inhomogeneous ones. The validity of this truncation scheme can be examined as a backreaction 
problem. One can view the homogeneous and inhomogeneous universes as the low and high lying modes 
of excitation of spacetime (this notion can be made more precise in terms of gravitational perturbations off 
a background spacetime as described by the Lifshitz wave operator), treat the former as constituting the 
system field (minisuperspace) and the latter as the environmental field. The magnitude of the backreaction 
of the latter on the former will give a measure of the validity of the minisuperspace approximation. This 
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approach was taken by Hu and Sinha who introduced the CTP CGEA |^5|. The same idea and method 
should be apphcable to backreaction effects of hard thermal loops on the soft sector in quark-gluon plasma 




In the high-low mode division, one can think of the system as the average of the field over a spatial volume 
A"'^, Ac = J being the inverse of the critical wavelength that separates the system from the environment. The 
Euclidean version of this coarse-graining leads to the so-called average effective action (AEA), introduced 
by Wetterich These two formalisms, Euclidean and real time CGEA, though sharing similar ideas have 
rather different ranges of applications. Wetterich's average effective action averages the fields over a space- 
time Euclidean region, while Hu's |^ coarse-grained effective action averages field over a spatial region. The 
Euclidean AEA is more suitable for near-equilibrium systems, and has indeed been applied to QCD problems 
like quark-gluon plasma in a quasi-stationary state, while the real-time CTP CGEA gives causal evolution 
equations for expectation values of quantum operator and is more suitable for the analysis of dynamical 
problems. We will invoke the AEA in Section VIII A. 



B. 'In-Out' Coarse Grained Effective Action 



We now introduce the m-oitt version of the CGEA. This is because is pedagogically simpler, and is closer to 
what one learns from textbooks, i.e., the Schwingcr-DeWitt effective action. Functionally it is also adequate 
for the introduction of ideas on coarse-graining and backreaction to the extent of seeing the running of 
coupling constants leading to a RG theory. But the effective equation of motion with backreaction obtained 
from an in-out version contains non real and non causal terms. This was first noticed by DeWitt and 
Jordan ||,|9| and Calzetta and Hu in a study of the backreaction of particle creation on a background 
spacetime or field. For a correct treatment of the backreaction of the environment on the system one needs 
the in-in or CTP version, which we will discuss in Sec. 6. 

Consider the model with action (III.l) in a spatially-flat Robertson- Walker universe. The scale factor 



a{t) = e" is left unspecified for now. Rather than viewing a(t) as a dynamical function determined from 
Einstein's equations we will, in the spirit of scaling as explained in an earlier section, regard a as a constant 
and t as a parameter. Although the space can have different scales a{t) at different times t, spacetime remains 
fiat at all times. (One can rescale x at different times to make them all equal in value.) Thus we can simply 
think of doing Minkowski field theory here, leaving a(t) as a parameter while carrying out coarse-graining 
in the 3-dimensional space. The content of this subsection can thus be used without reference to cosmology 
(simply set a{t) = 1 ). However we do want to tag the scale factor a{t) along so that later we can perform 
scaling transformations and discuss inflationary cosmology. 



For the RW spacetime, the action (III.l) reads 



S[(p] ^ J d^x J dtait)[-^(t>id2 + di - A + m^)(j) - ^A0^] (IV.l) 

where 

d2 = a^{t)d^/dt'^, di ^ a^{t)3Hd/dt, A = a'^{t)A (IV.2) 

are the second and first order time derivative operators and A is the spatial Helmholtz operator respectively. 
Here H = d/a — a and 

(t) = (t) (jr? -\- ^R) (IV.3) 



~X{t) = a2(i)A (IV.4) 

are the conformally-related generalized mass and field self-coupling parameter respectively. We will now 
assume di — d2 — 0, i.e., consider the case of eternal infiation. In this 'static' limit, one can introduce the 
momentum space representation on the spatial section 

0(f,<) = (2^)-3/2 [ d'^ke'^-^(l){k,t) (IV.5) 
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Then 

S[(t)] = J dta{t) j d^k[~^{k^ + m^{t))<j){k)<j){~k)] 

+ {2tt)-^ J dta{t) j d^kid^k2d^k3d^k4[-^X{t)d{ki + fca + ^3 + ki) 

4>{k , t)cp{k2 , t)cp{%, t)4>{ki, t)] (IV.6) 

For interacting fields in a static spacetime mode-mixing is manifest. If di 7^ "time" translation- invariance 
is lost and there will be frequency mixing. This becomes relevant in the context of critical dynamics, as 
mentioned in Sec. 2. 

Now we separate the field 4>{x, t) into two parts, 0< and (/)>, which one can refer to as the system and the 
environment, 



We assume that 0< contains the lower k wave modes and 0> the higher k modes and consider these two 
cases: 

Case A (critical phenomena) 

</)< :| k |< A/s, 0> : A/s <| k |< A, (IV.7) 

Case B (stochastic inflation) 

4)<_:\k\<eHa, :] k \> eHa, e w 1. (IV.8) 

Here A is the ultraviolet cutoff and s > 1 is the coarse-graining parameter which gives the fraction of total 
k modes counted in the environment. The separation of (j) can also be made in other manners, depending on 
the physical setup of the problem and the questions one asks. The formalism we present here can be easily 
adapted to other types of decomposition. 



S[(j)\ in ([V.6) can be written as 



S[(l)]^ J dta{t) J d^k{-^){P +m'^{t)(t)<{k,t)(l}<{~k,t)) 
+ J dta{t) j d^q{-]^){(f + m^{t))dpy{q,t)d^^{^q,t) 
-(-(27r)~3 I dta{t) I d^k[d^k2d^kid^k4-^X{t)) 



X(t)<^{ki,t)<j)<{k2,t)<j)<{k3,t)(t><{ki,t)S{ki + /c2 + ks + k^) 
+ {2t:)-^ J dta{t) j d^kid^k2d^k3d^q-{{-^X{t)) 

x4(t)^{ki,t)(j)<{k2,t)(t><:{k3,t)(t)y{qi,t)6{ki + h + h + qi) 
+ {2tt)-^ J dta{t) j d^kid^k2d^q-[d^q'2i-^Xit)) 

x6(?!)<(fci, i)0<(A:2, t)(/>>(gi, i)0>(g2, t)6{ki + k2 + qi + 9*2) 
+ {2Tr)-^ J dta{t) j d^kid\-{d''q-'2d^qli{-^X{t)) 

'X4:(f><{h,t)(p>{qi,t)(l)y{q2,t)(p>{q3)6{ki + qi + q2 + (fa) 
-t-(2^)-3 I dta{t) J d^q-id^q2d^q-'3d^q-l{-^Xit)) 

x0>(gi,i)0>(g2,t)0>(g3,t)0>('f4,i)<5(^L +92 +93 + 94) (IV.9) 
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where all integrals of k and q have the limits indicated in (IV.?) or (iV.S) depending on the cases in question. 
Henceforth we shall use the shorthand notation 



and denote 



G<(1) = [fc2 + rii^(t) - te]-\ G>(i) = [f + m^{t) - ie]-\ 



We can separate terms in ([V.9) into three groups 

= 5[0<]+5o[(/.>] + 5,[(/.<,0>] 

where 



(IV. 10) 
(IV. 11) 
(IV. 12) 



1 



1 r 



5['^<] = --G<i(l)0<(l)(/.<(-l) - -A(i)0<(l)0<(2)0<(3)(/.<(4),5(1 + 2 + 3 + 4) 



5o[0>] = -xG>i(i)'^>(i)'^>(-i) 



and 



5/[<^<, (A>] = -iA(i)[40<(l)0<(2)0<(3)0>(i)(5(l + 2 + 3 + 1) 

+6</.<(l)0<(2)0>(i)0>(2),5(l + 2 + 1 + 2) 
+40<(l)0>(i)0>(2)^>(3),5(l + 1 + 2 + 3) 
+0> (i)0> (2)0> (3)0> (4),5(i + 2 + 3 + 4)] 
= 5*1 + 5*2 + 5*3 + 54 



(IV.13) 



Here integration over k's and q's with the respective ranges is understood, and we have used the simplified 
notation of 1, 1 to denote fci, (fi, etc. It is easier to use Feynman diagrams for the propagators and vertices 
corresponding to the two fields (/>< and </>>. They are listed in Figure 2. 

The functional integral formalism for quantum fields is set up in the usual way. Assume that ai t — ±oo, 
the interaction is turned off so that the in-vacuum |0_) and the out-vacuum 10+) are defined. The vacuum 
persistence amplitude on the generating functional (0+|0_) gives the probability amplitude of the field 
remaining in the vacuum after the interacting system evolves in time. 

With the splitting = (/)<+ 0>, the Schwinger-DeWitt (in-out) generating functional becomes, 



Z[(f)] = N / L>0e 



,{jS[0<]+So[0>]+S/[0<,0>]} 



(IV. 14) 



Introducing an integral over of expiSo[0>] as the norm for the functional average in the denominator 
and changing the constant N to N' accordingly, we get 



(IV. 15) 



Here ( defines averaging over the </>> field under the free action 5o [</>>] (We write ( )o interchangeably 

to emphasize that the averaging is with the free action 5*0 [(?!)>].) 
Denoting 



Z[<j)] =N' J D(t)< e*^["*<l(y D(t)> e'So[0>]g»SH</><,0>]/ j i:»0>e'^"['^>l) 



(IV. 16) 



we get 
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Z[<I><]=N' J expire// 



The coarse-grained effective action is given by 

5e//[0<]=5[(/.<]+A5[0<], 

where 



(IV.17) 



(IV. 18) 



(IV. 19) 



If the interaction between the 0< and (/)> fields is smah (e.g. A ^ 1 in a At/)"* theory), one can expand this 
in a Dyson- Feynman series 



71=0 



(IV.20) 



where in the second hne terms up to the second order in Sj are written out exphcitly. 

Let us now examine each te rm in ([V.20) in detail, starting with the first order term < S"/ [</<<, (/)>] >0>. 
We see that terms in ([V.13) involve the average of products of (/>> fields, but those containing an odd 
numbe r of fields in the product average to zero. Also any product containing four (/)> fields (last term in 
( tV.13 )) averages out to a quantity which is independent of (/><, i.e., ((/)>(l))o = ((/)> (1)0> (2)0> (3))o = 
and (0> (!)(/()> (2)0> (3)0> (4))o is independent of (/)<(fc, t). Thus the only contribution to (S'/[0< ,</>>]) o is the 
quadratic product 



(0>(i)0>(2))o = -za-i(t)G>(i)<5(i + 2)S{h 



(IV.21) 



We denote by 5'/[0<, 4>>]j the terms containing j 0> fields, i.e. Si i (/)<(^<(^<(?!)>, S"/ 2 ~ 4><4><4>>4'>i etc. 
In the Feynman diagram depiction (Fig. 2c), the averaging over (/>> amounts to closing the lines representing 
the environment field propagators on themselves, and it is easy to see how the above conclusion is reached. 
Explicitly, the only contributing first order term is 

(5/[0<,(/)>]2)o = (27r)-3y j d^kid^g J d^qid^q26{-^)a{tyX{t) 

X(/)<(A;i,i)0<(fc2,<)((/)>((?'i,t)(/)>((?'2,i))o'5(fci + h + qi +92) 



Using ( |IV.2lD , we get 

(S"/ [(/)<, 0>]2)o 

where 

Sm^lt, s)i = — 



a(t)dt I d^k{--)(j)<{k,t) 5fh'^{t,s)i (j)<,{k,t) 



d^q [(f 



(IV.22) 



(IV.23) 



This is depicted as (1) in Fig. 2c. The subscript 1 on (5m^ denotes that it is the first order correction to m^. 
The second order correction is given by (5) an d (6) o f Fig. 2c and computed below. 

Now we examine the second order terms in ( [V.20 ). The second term in the square bracket (5'7[0<, 0>])q 
gives three non-contributing disconnected graphs. The first term ((5'7[0<, </'>])^)o has the following compo- 
nents: 

{Sho = {S!)o + {Sl)o + {Sl)o + {Sl)o 

+2[{SiS2)o + {8183)0 + {818^)0 
+ {8283)0 + {828^)0 + {838^)0] 



where we have written in short 



14 



S, ^ Si[<l><A>]3 0' = 1,2,3,4) 



as in Eq. (IV. 13) . In terms of the Feynman diagrams in Fig. 2c constructed from Fig. 2b, there are four 
contributing connected diagrams of the second order. They are denoted by (3)-(6) in Fig. 2c. We see that 
diagram (1) is the first order and (5) (6) are the second order corrections to the mass m^, (4) is a second 
order correction to the coupling constant A, while (3) generates (j)^ correction. 

Adding these contributions together we obtain finally for the coarse-grained effective action up to second 
order in A : 



Seff[^<\ = / a{t)dt / d^fc(-i)0<(fc,i)[fc2 + m2(t,s)]</)<(-fc,t) 

1 .r. 



a{t)dt J d^ki(Fk2(Fk3(Pki{~-)\{t,s,{k}) 

X0<(fci,i)0<(fc2,i)(/><(fc3,i)0<(fc4,i)5(fci + fc2 + 4 + ki) 
+ I a{t)dt I d^ki...d^k(i{-^)ve{t,s,{k}) 



6!' 



x0<(fci,i)...0<(fc6,i)(5(fci 



(IV.24) 



where the integration range of ki is as indicated in ( [V.7 ) and ( [V.8 ): | k \< A/s for the case of critical 

phenomena and | fc |< eHa for stochastic inflation. 

Here, 



m^{t, s) = rfi^(i) + Sm^{t, s)i + 5fh^{t, s)2a + 5rh^{t, s) 



2b 



A(t,s,{fc}) = ~Xit)+S~X{t,s,{k})i 



(IV.25) 



The subscripts under Srh^ and SX denote the order of correction to the mass and the coupling constants 
arising from averaging over fields. They are given by: 



1 



Sm'it,sha = {-::)mO)Y[X{t)/a{t)Y{27r)-' / d'' q,[cft + (t) - t 



(IV.26) 



Sm^t,s)2i, = (--)[i(5(0)]2[A(t)/a(t)]2(2^)-3 / d^q,i2n)-^ / d^q^i^f, + rh'it) ^ ze]-' 



[gf + m'^it) - ie]-^[{k - qi - q2f + m'^it) - it]' 



(IV.27) 



5X(t,s,{k}) = [«(5(0)] -a{t)-'X\t)(2T:y'' / d'^q'' + m\t) ~ le 



[{ki+k2-qf +m'^{t)+ief^ 



VQ{t,s,{k}) = -lQX^{t)[{ki+k2 + hfm'^{t) - it]-^ 



(IV.28) 
(IV.29) 
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V. BACKREACTION IN THE INFLATIONARY UNIVERSE: RENORMALIZATION GROUP 
EQUATIONS AND THE RUNNING OF COUPLING CONSTANTS 



We have so far discussed the separation and averaging processes leading to a coarse-grained effective action 
which contains the averaged effect of the environment on the system. Now we will introduce a rescaling of the 
k-space (truncated in critical phenomena, enlarged in stochastic inflation) to bring it back to the original size 
(but with less content in CP, more in SI). By requiring that this transformed Lagrangian has the same form 
as the one we started out with we obtain a set of renormalization group equations for the field parameters. 
At this point our earlier observation on the equivalence of scaling and inflation for the fleld and spacetime 
in question becomes relevant. That is, rescaling in a static picture serves the function of inflation in a 
dynamic picture. Our derivation here of the differential renormalization group equations is similar to the 
Wilson-Fisher or the Wegner-Houghton methods |Q . 

After introducing a rescaling of the k-space (spatial dimension D=3), 



k' = sk, (p'ik' , t) 



-^<P<ik,t) 



(V.l) 



the effective action (IV. 24) becomes 



SeffW] = / dta{t) / d^k'{~-)<l)'{k',t)[k'^ + s^m^{t,s)]<P'{k',t) 



j dta{t) J ( 



+ / dtait) I dk[ - ■ ■dk'^{~-)s^-°~X{t,s,{k'})(l)\k[,t)- ■ ■<P'{k'^,t)S{k[ + g + k'3^ 



+ I dta{t) I dk[ - ■ ■dk'^[--y-^^y^{t,s,{k'})c^\k[,t) - ■ ■cj,\k'^,t)6{k[ 



(V.2) 



Note that the fc'-integrations now resume the full range as k in ( IV. 7 ) or ( IV. 8 ). In terms of the rescaled 
variables (j)' , Scff will have the same form up to a certain order as the original Sog in terms of (p provided 
that we identify 



X'{t,s) 



4-D 



X{t,s) 



(V.3) 



The effective mass and coupling constants are given in (IV.25-IV.29). 

We can now proceed to tackle the coarse-grained integrals therein. At this point one needs to stipulate 
the system and bath separation such as the Cases A and B given as examples in Section 3. We will work out 
the details for critical phenomena (Case A) here. The case for stochastic inflation (Case B) can be obtained 
via the simple relation between Case A and B also mentioned in Section 3. For small changes in s ~ 1 -|- da, 
one can derive a set of differential renormalization group equations for the mass and the coupling constant 
as follows: 



dx/da — 2x — —y/{l + x) 



dy/da^ey+-y^/{l + xY 



(V.4) 



where 



X = 771^ (t, cr)A 



y = A(t, a)k-''na-^{t) 



(V.5) 
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Here, e = A — D, and is the solid angle integration in D-dimension. 

As expected, these equations have the same form as in the Ginzburg-Landau- Wilson model [|9| |o| || 
. They govern the flow of the field parameters (m^, A) as the scaling changes. There exist certain points in 
this parameter space known as the fixed points, where further application of the RG transformation leads 
to an invariant result. The fixed points are thus the steady-state solutions to the RG equations dx/da = 
and dy/da = 0. For the 0^ theory, it is well-known that there is a trivial fixed point at 



0, 



and a nontrivial fixed point at 







(V.6) 



= -e/6, y) = -2e/3 



Near the trivial fixed point the solution to (V.4) is 



(V.7) 



A 



„2<T 



B 



2(2 -e) 



(V.8) 



The critical point of interest to us is the nontrivial fixed point. We want to find out how in its neighborhood 
the field parameters flow towards this fixed point. By setting x = x*^ + Aa; and y = y*f + Ay we get 



d(Ax)/da = (2 



1 



1 



1 



and 



d{Ay)/da = -eAy 



(V.9) 



with the solution 



Ax\ 
Ay j 



= A 



= (2-e/3)(T 



B 



4(1 + e/6) 



(V.IO) 



A. Cosmological Consequences 



As we explained earlier, in eternal inflation the renormalization group transformation parameterized by 
s is equivalent to time evolution parametrized by t. Thus under such conditions we can treat the coupling 
constants like X(t, s) as depending on either s (flow) or t (evolution). The fixed point corresponds to the 
extreme infrared limit which in the case of eternal inflation lies at infinitely late time tf. In practice it 
corresponds to the time the inflation regime ends ti ~ tf before reheating sets in. Thus in (V.IO) Ay is a 
measure of how y changes from t to ti. Explicitly, with e = 1 and a = J Hdt, the 0(e) expansion result 
gives 



(V.ll) 



At the fixed point tf, y{tf) = -|e. Recall that y = \{t)A^''na^^{t) and X{t) = X{t)a'^{t) ,this yields 



l-7Be"^'-'^^ 



(V.12) 



Comparing X{t) at two times t > to, 



Xito) 



-H(t^to) 



1 - 7B{e 



(V.13) 



Assuming tf >> t,to (which is clearly satisfied for realistic inflation tf > toe^^) the second term in the 
square bracket can be neglected. We get. 
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\{t)/\{to)~e-"'-'-'°^ (V.14) 

Thus the couphng constant A of the effective field theory (in a system which is depleting in content) actually 
decreases with time. |^ This is a new effect which could have important consequence for the galaxy formation 
problem. 

As we know from standard calculations the density contrast Sp/p depends strongly on A. For GUT 
processes one needs to assume an unnaturally weak A ~ 10^^^ in order that Sp/p « 10~^ at the time the 



fluctuations reenter the (RW) horizon. What (V.14) implies is that the strength of A at time Iq when the 
fluctuation mode corresponding to galaxy scale left the de Sitter horizon is much weaker than X a.t f^, the 
time the Hubble size left the de Sitter horizon. The running of the coupling constants seems to provide a 
way to reduce their strength. The reason behind this phenomena is, as we recall, due to the backreaction of 
the environment of high k modes on the system of lower modes, which in this particular choice is diminishing 
in time. 

In contrast, for stochastic inflation where the system (consisting of modes with physical wavelength greater 
than the de Sitter horizon) increases in content, the coupling constants will, from this analysis, actually 
increase exponentially with time during the inflationary regime (era from to to ti). Because a formal 
analogy exists between this case (Case B) of stochastic inflation (SI) and critical phenomena (CP, Case A) 
one can indeed just read off the result from the results we obtained earlier. The transcription is simply i) 
Changing A in CP to eH in SI. Recall that A refers to the ultraviolet cutoff in CP, whereas is the 

horizon size. The ultraviolet cutoff for SI is assumed to be oo. ii) Changing s in CP to s^^ in SI. Recall 
that the coarse-graining parameter s plays the role of the scale function a{t). Thus in CP is now acting 



like a{t) in SI [see Eq. ( IV.S )]. Note that s > 1 in CP but s < 1 in SI. Note also that transformation (ii) 
is equiva lent t o a time reversal, i.e., changing At —At. The way A runs in stochastic inflation is given 
by Eq. ( |V.12 ), with —t replacing t there. The ensuing discussion of the implications on galaxy formation 



is similar, except that, of course, one reaches the opposite conclusion, i.e., the strengthening of A in time 
increases the density contrast as each mode reenters the RW horizon. 

The opposite results arising from these two cases should not be viewed as unsettling. They arise from very 
different stipulations of the system (CP decreasing in time, SI increasing). The cosmological consequences 
will depend on whether and how A runs, but the theoretical issues this investigation brings out should be 
valid. What is shown here is that if one takes into account the interaction between the environment and 
the system, the running of coupling parameters in the theory as observed in the system is an unavoidable 
consequence. How they run (with time or scale or energy) will depend sensitively on how the system is 
selected with respect to the environment, and, of course, the nature of the interaction (gauge fields would 
presumably run differently from X<p'^, as their UV and IR behavior differ). 

On this point, let us reexamine the physical criteria for the choice of the system versus the environment in 
our examples. The cutoff is determined by the highest k mode which left the de Sitter horizon, which in turn 
depends on the duration of inflation. At the start of inflation of course there is no knowledge when inflation 
will end and the modes do not know which ones among them are to be counted in (the system) and which 
are to be out (the environment). A more straightforward division would be to have a flxed cutoff at the 
beginning of inflation, e.g., those with k greater than the galaxy size ( ka ) constitute the environment and 
those smaller than kc but greater than be the system. This would give rise to running, but presumably 
at a different rate. Indeed, if one is interested in the behavior of a particular mode, say, that of the galaxy 
scale kc, one can regard this one mode as the system and the rest as the environment. At all times this 
one mode is interacting with the rest of the whole spectrum, some of these leave the DS horizon earlier, 
some later. The behavior of this particular mode determines at a critical time t^ (when it leaves the DS 
horizon) its own amplitude when it reenters the RW horizon, but it is influenced by all the lower and higher 
modes at earlier times. The effect of coarse-graining for this choice of the system (single mode) is similar to 
the minisuperspace problem in quantum cosmology. It is expected that the coupling constants will not run 
in these circumstances. Even in more general cases it may also be that the changing A as experienced by 



^Of course, even though the running coupling goes like A(t)/A(to) — e"^'*"*"', a loop expansion of a quantity such 
as the density- density correlation function will show that in the same limit each loop appears with a factor of e^'*"*"' . 
The self-consistent cancellation of these two factors is associated with a crossover phenomenon and in principle could 
be accounted for by using an H-dependent (environmental-friendly) RG [^. We thank C. R. Stephens for this 
remark. 
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modes in the system may not exert any effect on galaxy formation, because at the moment any kth mode 
departs the de Sitter horizon it will have the same amplitude no matter what history of interaction it has 
experienced, and it is this amplitude which determines the density contrast for galaxy formation in the RW 
era. This point, however, raises an interesting theoretical issue. 

B. Theoretical Implications 

We see that different choices of the integration range in the calculation of the CGEA give rise to very 
different behavior of the coupling constants. In critical phenomena we stipulated the system to be in the 
same way as the real space renormalization group treatment of critical phenomena, i.e., that the ultraviolet 
cutoff is fixed at the start of inflation- this corresponds to the inverse of the lattice size, which is a fixed 
value for condensed matter considerations. In each iteration the UV cutoff is redshifted, making the effective 
range of frequency smaller at successive periods. In critical phenomena scaling is an artificial procedure to 
facilitate the approach to the infrared limit. However, in inflation every different time interval corresponds 
to a realistic physical situation. If one adopts the coarse-graining scheme similar to critical phenomena, one 
would get a diminishing range for the system as the universe expands. This is the cause of the changing 
coupling constant, and is also the source of a serious problem. Imagine that if one turns the problem around 
by demanding that the physics observed in each of the subsequent moments should be identical, specifically, 
if it possesses the full range from the Planck scale to infinity (k from to inverse Planck length)- this is 
certainly what one usually assumes, otherwise the physical range we observe in the post-inflationary era 
such as in today's universe would only consists of the low k modes-if one makes this reasonable assertion, 
then we would not have the running coupling constants problem. This is equivalent to assuming that the 
physical wavemodes p = k/a rather than the intrinsic wavemodes k has the integration range between zero 
and the inverse Planck length. However, making such an adjustment in order to make the range identical at 
each subsequent moment would require a mechanism to 'replenish' the high frequency modes between the 
redshifted UV limit and the Planck length, which would open up another serious problem. Notice that this 
problem is not particular to inflation nor curved spacetime. It is already there for any theory in a dynamical 
setting, e.g., in cosmology. It only becomes more serious in inflation because the redshifting is exponential. 
Similar questions would arise in black holes when one tries to compare what is observed between different 
observers from the event horizon to asymptotic infinity. The relation between black hole and inflation can 
be easily understood by viewing them in terms of exponential scaling, as in the dynamical finite size effects 

One can analyze this problem in the context of the renormalization group theory by using an explicit 
cutoff (defined at one instant as the Planck scale) . By demanding that this cutoff be the same at all instants 
and at all energies, one would presumably not get any flow. This would in turn demand some theory on the 
physical consequence of transPlanckian frequencies 

PART THREE. Backreaction and Noise from CTP- CGEA 

VI. THE CLOSED-TIME-PATH COARSE- GRAINED EFFECTIVE ACTION 

As we introduced the CGEA in the 'in-out' version Q] we remarked that for a correct treatment of 
backreaction problems, the 'in-in' or closed-time-path CTP version of CGEA is the right way to proceed, as it 
gives a real and causal equation of motion for the effective dynamics of the open system. The CTP CGEA was 
first introduced by Hu and Sinha ||35|| to analyze the validity of the minisuperspace approximation in quantum 
cosmology. The influence functional method [E5[ for interacting quantum fields (based on quantum Brownian 



^This issue raised originally in the context of inflationary cosmology in |Q has a parallel in black holes, which was 
pursued independently by Jacobson, Unruh and others. For a discussion of trans-Planckian modes problems and 
explorations for its implication for quantum gravity, see, e.g., H. 
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motion |Q) was introduced by Hu, Paz and Zhang in and Calzetta and Hu ||^, and applied to a 

range of problems in semiclassical gravity and cosmology such as decoherence, structure formation, entropy 
generation and reheating p^ , p3| , p8 
CGEA for a self-interacting theory in a Robertson- Walker spacetime to discuss the quantum to classical 
transition of the low modes of the field, and to put on firmer grounds the quantum theory of structure 
formation in inflationary models. Greiner and Miiller also used the CTP CGEA to analyze the classical 
limit of the soft modes of a quantum field when the hard modes are in thermal equilibrium. Dalvit and 
Mazzitelli ||5^ showed that it is possible to write an exact equation for the dependence of the CGEA on the 
scale that separates the system and the environment. This is expected to be useful for a nonperturbative 
calculation of the CGEA, and to discuss the appearance of noise in the renormalization group equations. 

In this part we will review the above mentioned works. We will start with the necessary definitions and 
the perturbative evaluation of the CTP CGEA for a Xcj)'^ scalar field in a Friedmann-Robertson- Walker 
(FRW) background spacetime. We will describe how to use the CTP CGEA to analyze the validity of the 
minisuperspace approximation in quantum cosmology and the issue of quantum to classical transition of 
fluctuations which give rise to structure formation in inflationary models. We will finally describe attempts 
to compute the CTP CGEA using non perturbative approximations. 



49 1 . Lombardo and MazziteUi Q computed perturbatively the CTP 



A. Perturbative evaluation of the CTP CGEA 



We consider again the scalar field action given by Eq.(III.l). In a fiat Robertson- Walker spacetime with 
metric ds^ — dt^ — a'^{t)dx^ = a?'{rf)[drf — dx^], the action can be written as 

S{a, x) = 1 d^x[\ri^^d,xd.X - - \)Ra\^ " ^x'] (VLl) 

where x = o,(t>. From now on, d^a; denotes d^x drj. To see the fiat space results simply set a = 1 and i? = 0. 
Let us make a system-environment field splitting 

x{x)^X<{x)+x>{x), (VI.2) 

where we define the system by 

f d^k ^ 

X<{x,t) ^ -——^(t>{k,rj) exp ik.x, (VI.3) 

and the environment by 

[ d^k 

X>{x,v)^ -— r(/)(/c,77)expi/c.f. (VI.4) 

J\k\>A^ (271-)^ 

The system- field contains the modes with wavelengths longer than the critical value A^^, while the bath or 
environment-field contains wavelengths shorter than A~^. Ac corresponds to of previous sections. 
After the splitting, the total action ( VI.l[ ) can be written as 



S[a,x] = Sq[x<] + Sq[x>] + 5'„t[a,x<,X>], (VI.5) 
where So denotes the kinetic term and the interaction part is given by 

f A fM'^ 2 ^ 4 2 A 4 

Sint[a,x<,X>\ = ^ ^i—X< + ■^X< + —X> + + 

\x\{-)xl + ^x\x> + lx<xl}- (VI.6) 

with IVP = m^a? -f - \)Ra^. 

We are interested in the influence of the environment on the evolution of the system. Therefore the CTP 
CGEA S^[a,x<^o! ,x'<\ is the object of relevance. It is defined by 
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r rx>s r-x>f 

expj5'Aja,X<,a',x<] = expi(5o[x<] - -Soixy) / dx>f j Vx> / 2?X> 

X expi{S'o[x>] + S^nt[a,x<,X>] - Sq[x'>] - S^nt[a' ,x'<,x'>]}- 



(VI. 7) 



The integration is over all fields x> (and x>) with positive (and negative) frequency modes in the remote 
past that coincide at the final time x> = X> = X> f 

We will calculate the CTP CGEA perturbatively in A and M^, up to quadratic order in both quantities. 
After a simple calculation we obtain 

'5'Aja,X<,a',X<] = 5'o[x<] - 5'o[x<] + {{Stnt[a,X<,X>])o - {Sint[a' ,x'<,x'>])o} 
+ '^{{Sfnt["',X<,X>])o - [('5'j„t[a,X<,X>]>o]^} 

-^{{S^nt[a, X<,X>]Stnt[a, x<,X>])o - {St„t[a, x<,X>])o{Sz„t[a , Xo X>])o} 

+ '^{{SLW,x'<,x'>])o-[{S.nt[a'x'<,x'>U'}, (VI.8) 

where the quantum average of a functional of the fields is defined with respect to the free action 5*0 

{B[x>,X>])o= J dx>fp"'vx> Vx'yexp{So[x>]~ So[x>]}B. (VI.9) 



Equation (VI.8) is the in-in version of the Dyson-Feynman series 
We define the propagators of the environment field as 

{x>{x)x>{y))o = iG^+{x-y), 



(VI. 10) 



{x>ix)x'>{y))o^-^GlLix-y), 



(vi.ii) 



{x'>i^)x'>iy))o = -^G'^-ix-y)■ 



(VI. 12) 



These propagators are not the usual Feynman, positive-frequency Wightman, and Dyson propagators of the 
scalar field since, in this case, the momentum integration is restricted by the presence of the (infrared) cutoff 
Ac. The explicit expressions are 



G+^{x~y) 



|p|>A< 



p'' + ie 



(VI. 13) 



G^^_(x-y) = - 



Ip1>Ac 



d-^p 
(2^ 



e*f(^-^)27^^(5(p2)e(/), 



(VI. 14) 



G^_(x-y) 



|P|>A< 



d'^p 

(2^' 



p^ — ie 



(VI. 15) 



As an example, we show the expression for the propagator G^!^. The usual Feynman propagator is 

G++(x) = -\l-^<5(a), (VI.16) 



while 



GX\{x) — 3-3- 



cos[Kc{r - x°)\ cos[Kc{r + x^)] 



r[r — x") 



r[r + a;0) 



1 y sin[Kc{r — x'^)] sin[Kc{r + x^)] 
8^ 



r(r — a;0) 



r(r + a;0) 



G++{x)-Gf<'''{x) 



(VI. 17) 
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where a — ^x"^ 



The CTP CGEA can be computed from Eqs. (|VL^)-( |VI.l^ ) using standard techniques. After some 
algebra we find 



SAAa,x<,a',x'<] = So{x<) - So{x'<) - ^ / d^x{M^x)x% - M 



^(X<(^) - XM) + -iGl^mU^ix) ~ M\x)) 



d^x / d^y 



^x\{x)Gl^+{x - y)x\{y) - ^x\{x)G^'^{x - y)xl{y) 



-^x'i(^)G^_(x - y)x'<^(y) - \M\x)zGli{x - y)M\y) 



+ ]^M\x)iGX'^{x - v)M\y) - -^M\x)iG'^^^{x - y)M\v) 

+^x<(^)G^+'(x - y)x<{y) + ^x<{x)G^+-{x - y)x'<iy) 
X 



^gX'<(^)G^^l'(x-2/)x'<(2/) 

where we introduced the notation = AP 
Defining 



P± = i(xUx'<') 



R± = -{x%±x^) 



(VI. 18) 



x± = ^{x<±x'<) 



and using simple identities for the propagators, the real and imaginary parts of the CTP CGEA can be 
written as 

ReSA^ = 5o(x<) - So{x'<) - Ay d^x{^P^{x) + ^G^=+(0)Q_(x)} 

J d^x J d^yQ{y''^x°){-^R+{x)ReGlx{x-y)R-{y) + 

^Q+{x)ImG'l'^ix ~ y)Q^{y) + ^-x+{x)ReGXi{x - y)x-{y)}. (VI.19) 
ImSA^^X^ J d^x J dM-^RAx)lTT^Gl^_^{x ~ y)R_(y) 

-^Q^ix)ReGl'^ix - y)Q^iy) + ^x-{x)ImGi^^ix - y)x-(y)}. (VI.20) 

The real part of th e CTP CGEA in Eq. ( |VI.19| ) contains diver gences and must be renormalized. As the 
propagators in Eqs. ( VI.10)-( VI.12) differ from the usual ones only by the presence of the infrared cutoff, 
the ultraviolet divergences coincide with those of the usual Ax^-theory. The effective action can therefore be 
renormalized using the standard procedure. 

Consider the square of the Feynman propagator. Using dimensional regularization we find 



Gl^^ix) = Gl+ix) + G^T'^^'^i^) 2G++{x)G'_^--^'ix), 



.(b1<Ae)2 



(IpKAc), 



where 



GlAx) 



1 



+ -0(1) - 47r](5*(a;) + - ri{x) - Log[ATrfj,'^ 



(VI.21) 
(VI.22) 



j:{x) 



dSe'P^Log\p% 



(VL23) 
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''^''^ ^ (2^ i dV^e(p2). (VI.24) 
Note that the divergence is the usual one, i.e., proportional to —^S'^{x ~ y) and independent of Ac. Conse- 



quently, the term ImG^\{x — y)Q+{x)Q-{y) in Eq. ( VI.19| ) is also divergent and rcnormalizes the coupling 



constant A and the constants that appear in the gravitational action (as usual, in order to renormalize the 
theory of a quantum field in a curved space, it is necessary to include in the gravitational action the Einstein 
Hilbert term, a cosmological constant, and terms quadratic in the curvature tensor). The other divergences 
can be treated in a similar way. One can also check that the imaginary part of the effective action does not 
contain divergences]^ 

As we will show with the help of QBM models, the (nonlocal) real and imaginary parts of Sa^ [a, X< j ^'i X<] 
can be associated with the dissipation and noise respectively, and can be related by an integral equation 
known as the fluctuation-dissipation relation. 

Similar perturbative results for the CTP CGEA have been obtained by Greiner and Miiller in Ref. [Q. 
In order to derive effective field equations for the soft modes of the scalar field, they computed the CTP 
CGEA assuming that the modes in the environment are at thermal equilibrium at a temperature T such 



that Ac <C T. Their result for the effective action is essentially given by expressions (VL19) and ( VI.2C ), for 
the particular case a{t) = 1 and replacing the vacuum propagators by thermal propagators in order to take 
into account the state of the environment. 

In Section VII. B we will see that the field equations derived from the CTP CGEA are real and causal. 



B. Applications: Dynamics of System Fields/Modes with Backreaction of Environment Fields/Modes 

The coarse-grained effective action is a very useful method to treat coarse-graining and backreaction 
problems. Examples include stochastic inflation [H,^2|j5^ and reheating in inflationary cosmology p3| , ^ , 
effect of hard thermal loops in QCD plasma, and the effect of individual atoms on a BEG condensate. The 
first instances the CTP effective action was applied to were for backreaction of quantum fluctuations and 
particle creation on the background spacetime or for interacting quantum fields . The CTP CGEA 
was first introduced | |35| ] to address the validity of the so-called minisuperspace approximation. Although 
it was introduced in the context of quantum cosmology, the method has a wide range of applications. The 
relevant issues are the backreaction effect of the inhomogeneous modes in an interacting quantum field 
on the homogeneous mode; and the validity of quantizing a truncated theory: does it preserve the salient 
features of a fully quantized theory? These issues underlie many problems in physics, esp ecially when we 
view effective theories as playing a more fundamental role in the description of Nature |55|] . 



VII. MASTER AND LANGEVIN EQUATIONS IN QUANTUM FIELD THEORY 

We now show how to use the open system concepts and techniques to give a first principles derivation of 
the evolution equation for the reduced density matrix that describes the system under the influence of the 
environment. We will use the influence functional method to extract the noise and dissipation kernels. The 
open system framework is best illustrated by the quantum mechanical problem of Brownian motion (QBM). 
Feynman and Vernon flrst treated this problem with the influence functional method. Using this method, 
the master equation was derived by Caldeira and Leggett p6[ | for Markovian processes (Ohmic environment 
at high temperature), and by Hu Paz and Zhang |56| for a general environment including nonMarkovian 
processes. (See also and [Q). As we will soon find out, the influence functional method is intimately 
related to the CTP coarse-grained effective action [^9|j7|] . 



Of course, a successful ultraviolet renormalization does not guarantee that an approximation scheme such as RG- 
improved perturbation theory will be well behaved. A good example is in Section V, where loops depend on a factor 
gif(t-to) which would invalidate perturbation theory. Further 'infrared' //-dependent (environmentally friendly) 
renormalization of A is needed. 
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Hu, Paz and Zhang p7|,|6[ first generalized the quantum mechanical problem of Brownian motion to 
quantum fields, taking the system and environment as two independent fields. Lombardo and Mazzitelli ||5Ct| 
treated the same problem, taking as system and environment the low and high frequency modes of a single, 
self-interacting scalar field. We will follow their treatment here. 

We begin with a brief review of the problem of quantum Brownian motion. Denote by x the coordinate 
of the Brownian particle, by its frequency, and by qi the coordinates of the oscillators in the environment. 
The influence of the environment on the Brownian particle can be described by the reduced density matrix 
Pr{x,x',t) that is obtained from the full density matrix by integrating out the environmental degrees of 
freedom Qi. 

For a linear coupling xqi, the master equation for the reduced density matrix Pr{x, x' , t) is of the form | |6C[ | 
ihdtPrix,x',t) = {x\[H,pr]\x') - ij{t){x - x'){dx - dx')pr{x,x',t) 

+ f{t){x - X'){d^ + d^,)pr{x, X',t) - iD{t){x - X'fprix, x' , t) , (VII.l) 

where the coefficients j{t),D(t) and f{t) depend on the properties of the environment (temperature f3~^ 



and spectral density I{uj))- The first term on the RHS of Eq. (VII.l) gives the usual Liouvillian evolution. 
The second is a dissipative term with a time dependent dissipative coefficient j(t). The last two are diffusive 
terms. The one proportional to the anomalous diffusion coefficient f{t) does not produce decoherence, i.e., 
the off-diagonal terms in the density matrix are not suppressed in time. To see this in a simple example, 
assume that pr — poe^^'^*''^^^^ ^ . Inserting this into Eq. ( VII.l[ ) it is easy to prove that the /—term 



produces an oscillating function B{t) but no damping. (We refer the reader to Refs. [|60|,|6l| for a more 
detailed justification). The last term with diffusion coefficient D{t) 



D{t) = I ds cos{ns) I dw I{uj)coth{-l3huj) cos{ujs), (VII.2) 







gives the main contribution to decoherence. Indeed, an approximate solution of Eq. (VII.l) is 



Pr[x , x' ; t] ~ pr[x , x' , 0] exp — [x — x')"^ I D{s)ds , (VII. 3) 



and we see that the off-diagonal terms of the density matrix are suppressed as long as D{s)ds is large 
enough. For non- linear couplings like x"q^, one expects the master equation to contain terms of the form 
ii5("''")(t)(a;" -a;'")Vr- 

We now proceed to quantum field theory. The total density matrix (for the system and environment fields) 
is defined by 

P[X<,X>,X<,X>,V] = {x<X>\p\x<X>)i (VII.4) 

where |x<) and |x>) are the eigenstates of the field operators x< and x>, respectively. For simplicity, we 
will assume that the interaction is turned on at the initial time ijo and that, at this time, the system and 
the environment are not correlated. (The physical consequences of such a choice is elaborated in . More 
general initial conditions can be introduced by a stipulated preparation function |6^.) As such, the total 
density matrix can be written as the product of the density matrix operator for the system and for the bath 

p[Vo] = P< [Vo]p> [Vo]- (VII.5) 

We will further assume that the initial state of the environment is the vacuum. 
The reduced density matrix is defined by 



Pred[x<,X<,v] = J 1^X>P[X<,X>,X<,X>,V]- (VII.6) 
and evolves in time according to 

Pr[x</,X</,'7] = y" dx<i J dx<,Jr[x<f,X<f,il\x<t,X<i,'no]Pr[x<iX<i,'no], (VII.7) 
where Jr[t, 770] is the reduced evolution operator 
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Mx<f,x'<f,v\x<^,x'<^,VQ]= f"' 'Dxk T^x'<ewUs[x<]- S[x'<]}F[x<,x'<]- (Vlli 

The influence functional (or Feynman- Vernon functional) F[x<:,x'<] is defined as 

//■ r rx>f rx>f 

dx>i / dx'>,p>[x>^,x'>i,■no\ / dx>f / Vx> / X>x> 



X exp J{^[X>] + ^mt[x<, X>] - ^[X>] - ^mt[x<, X>]}- (Vn.9) 

We see that when the environment is initially in its vacuum state, the influence functional coincides with 
the CTP CGEAjllQ. Treatment of two self-interacting fields via the influence functional leading to noise 
is contained in |8| 

As the propagator Jj. is defined through a path integral, we can obtain an estimation using the saddle 
point approximation: 

J'r[x</,X</,?7lx<i,X<»,?7o] ~ exp^S'Ajx<,X<'], (VII.IO) 

where x<(x<0 is the solution of the equation of motion ^^^^1^^=^'^ = with boundary conditions 

X<('?o) = X<i(x<i) E^nd x<{v) — X</(x</)- This formula enables us to analyze the quantum to classical 
transition of the system field using the perturbative evaluation of the CTP CGEA described in the previous 
section. 



A. Master equation and decoherence of long wavelengths 

The decoherence effects are contained in the imaginary part of the CTP CGEA, which is already of 
order A^. As a consequence, in the evaluation of ImS\^ we can approximate x'< by the solution of the 
free field equation satisfying the appropriate boundary conditions. For simplicity we will consider massless, 
conformally coupled fields. Therefore the classical solution reads 

^, sinjkos) sm[fco(?7~s)] ^ ^, ^ 

^<(^'^) = [^<^im^ +^<^^m(Mr]^"^(^°-^^ = X<is)cos{ko.x), (VII.U) 

where we assumed that the system-field contains only one Fourier mode with k = k^. This is a sort of 
"minisuperspace" approximation for the system-field that will greatly simplify the calculations. 

As in the QBM problem, we can derive the diffusion coefficients from the master equation, which would 
provide relevant information about decoherence. Hu Paz and Zhang |^,^ derived the noise kernel of two 
interacting 0** quantum fields in de Sitter space and analyzed its behavior in relation to decoherence. Lom- 
bardo and Mazzitelli pO| did the same, with one self-interacting field split into two mode sectors. They 



followed the method proposed in |60|, and computed the time derivative of the propagator J^, eliminating 
the dependence on the initial field configurations x<i ^nd x'<i that enters through x'< snd x< ■ The master 
equation is of the form 

r - 1 (x^ ~ x''^ 

ina,,p^[x</,X<y,77] = (X</| Hren,Pr lx</)-^A^[ ^^^^52'^ Di{ko;r]) 

+ 32 D2{ko;T]) ^ D3{ko;r))]pr[x<f,X<f,v] + (VII.12) 

Due to the complexity of the equation, we only show the correction to the usual unitary evolution term 
coming from the noise kernels. 

This equation contains three time-dependent diffusion coefficients Di{'q). Up to one loop, only Di and D2 
survive and are given by 
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Di{ko;ri) ~ / ds cos^ (A:os)/mG^!^ (3fco ; ry — s) 

^0 



6fc, 



Jo 



ds cos^{kQs) cos(3fcos) 0(3fco — 



2fco'7 + 3sin{2kori) + |sm(4fco'7) + |sm(6A:o'7) 



A, 



576fc^ 



< /co < Ac 



(VII.13) 



D2{ko;r])^ / cos2(fcos)(i?eG+!,_^(2fco; 7? - s) + 2i?eG+f (0; 7? - s)). (VII.14) 
Jo 

Using that 

ReG+'I{2ko; jj-s)^ — I dp dzcosUp + z)s] 

ko L 7a, Ja, 

+ / dp/ dzcos[(p + z)s] L (VII.15) 

J2fco+Ac Jp-2ko ) 



ReGlf{0;f]-s) = 27t5{s) ~ 2^^^^^^^ k (VII.16) 



the D2 diffusion coefhcient can be written as 

D2{ko;v) = J - - ^)^*[2rKA, ~ fco)] 



-(2 - ^^)S^[2Acv] + ^)5*[2^(Ac + ^o)] - (f + ^)5z[2ry(2fco + A,)] 



cos 



[2Ac77] cos[2'i]{Ac + ko)] cos[2ry(Ac - fco)] cos[2?7(2fco + A^)] 



4fcot 4fco?7 4fco?7 4fco77 

where Si\z] de notes the sine-integral function 



(VII. f 7) 



Eq. (VII. 12) is the field-theoretical version of the QBM master equation we were looking for. In our case, 



the system is coupled in a nonlinear form. Owing to the existence of three interaction terms (x<X>7 X<X>j 
and X<X>) there are three diffusion coefficients in the master equation. The form of the coefficients is fixed 
by these couplings and by the particular choice of the quantum state of the environment. 



Our results are valid in the single-mode approximation of Eq. (VII.ll). In this approximation one obtains 
a reduced density matrix for each mode fco, and neglects the interaction between different system-modes. 
Due to this interaction, in the general case, pr will be different from Yikf, Pr{ka)- This point deserves further 
study. 

A detailed analysis of the quantum to classical transition in the model we are considering is a very 
complicated task. One should analyze in detail the master equation and see whether the off-diagonal elements 
of the reduced density matrix are suppressed or not. One should also study the form of the Wigner function, 
and see whether it predicts classical correlations or not |^,^ Having in mind the analogy with the QBM, 
here we will only be concerned with the diffusive terms of the master equation. By examining how the value 
of the diffusion coefhcients change in time, one can get an indication of how effective decoherence is [|6|,p0[. 
This simplified method can only provide a rough approximation. 

A mode in the system will decohere if the diffusion coefficients are different from zero during an appreciable 
period of time. Therefore we ask the following question: which is the maximum value of the cutoff Ac such 
that, a few e-foldings after the initial time, all modes with fcp < Ac still suffer the diffusive effects? The value 
of the diffusion coefficients at a given time depends on the value of the adimensional quantity I = Ac?7. For 
the particular case of a de Sitter spacetime we have a(t) — exp(_fft) and 

/ = A.„ = ^(l-i, 
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which, a few e-foldings after the initial time, is approximately given by I ~ Ac/ H . We can distinguish two 
possibilities. When I ^ 1 both coefficients Di and D2 are appreciably different from zero for all values of ko. 
On the other hand, when I ^ 1 the situation is completely different: Di is very small everywhere while D2 
is also very small in the infrared sector. The diffusive effects are important only for ko ^ Ac. The conclusion 
is that, in order to have decoherence for all modes with ko < Ac, we can include in our system only those 
modes with wavelength larger than H~^. This is consistent with Starobinsky's original suggestion. If we 
include wavelengths shorter than -ff^^, the frequency threshold in the environment increases, the infrared 
sector of the system cannot excite the environment and therefore does not decohere. 

A more realistic calculation should include a time dependent cutoff since the system should contain, 
at each time, the modes with kph — — < H . More importantly, the scalar field should be minimally coupled 

= 0) to the curvature. In spite of this, we think that our example illustrates the main aspects of the 
problem. Indeed, one could repeat the calculations for a g iven mode of the minimally coupled scalar field, 
and obtain a master equation, similar to our Eq. (VII. 12). The main difference would be that, for ^ = 0, 
the propagators in De Sitter spacetime would differ from their flat-spacetime counteparts Eqs. (VI. 13- VI. 15) 
(See Zhang's thesis [0). The evaluation of the diffusion coefficients using curved spacetime propagators 
would be more realistic for discussions of structure formation in inflationary cosmology. 



B. The Langevin equation 

We now show how to derive a stochastic (Langevin) equation for the system field from the CTP CGEA. This 
equation takes into account the three fundamental effects of the environment on the system: renormalization, 
dissipation and noise. 

The real and imaginary parts of the CTP CGEA for our model are given in Eqs. ( IVI.19D and ( |VI.20| ). 
One can regard the imaginary part of S'a^ as coming from three noise sources v{x), f (x), and rj{x) with a 
Gaussian functional probability distribution given by 



P[u{x),ax),T^{x)]=N,N^Nne,^v{ - ^ I j d^yv{x) 



-ImGl 



X exp 
X exp 



d'^x / d^yS^ix) 



d X d yrj{x) 



'2A ^A^3 



-1 



(VII. 18) 



where N^, iVj, and TV^ are normalization factors. Indeed, we can write the imaginary part of the influence 
action as three functional integrals over the Gaussian fields v{x), ^i^): ^^nd ri{x): 

Vv{x) j Vi{x) j I?77(x)P[i/,e,r7]exp-j|i?-(a;)i^(x)+g_(x)C(a;)+X-(a:)»7(a;) 



exp 



d^x / d^y 
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R-{x)ImG'lAx,y)R^{y) 



+^Q_{x)ReGil{x,y)Q^iy) - yX-(a^)/mGf+(x,y)x-(y) 



Therefore, the CTP-CGEA can be rewritten as 
where 



Seff [x< ,x<,>^,£„v]= RsSa, [x< , X<] 



d^x 



R^{x)v{x) + Q-{x)i{x) + x-{x)ii{x) 



(VII.19) 



i/n j VvP[v] j ViP[i] j 2?,7P[,7]exp|i5c//[x<,X<,^^,e,^]|, (VII.20) 



(VII.21) 



From this effective action it is easy to derive the stochastic field equation for the system 



27 



^MX<^^|^^^^_0 . (VII.22) 



It is given by 

□x< - 



+ ^X\<{x) J d^ye{xo-yo)ImGl'+{x,y)xl{y) + l>^^ J d^ye{x^ - y^)ReG^i{x,y)x<{y) 

= \v{x)xl{x) + £,{x)x<{x) + \^{x) . (VII.23) 

The field equation is real and causal, as expected. We see that it contains multiplicative and additive colored 
noise. The non linear coupling between modes makes the Langevin equation much more complicated than 
the usually assumed white noise equation (see Section III C). 

Greiner and Miiller obtained a similar stochastic equation in flat spacetime, for a thermal environment. 
They analyzed in detail the dissipative terms in the Langevin equation. In particular they found explicit 
expressions for momentum dependent dissipation coefficients using a Markovian approximation for the soft 
modes. 



VIII. RENORMALIZATION GROUP FROM CTP CGEA 



A. Towards a nonperturbative evaluation of the CTP CGEA: The exact RG equation 

The main drawback of the results presented in the previous sections is that we have evaluated the CTP 
CGEA only perturbatively. As is well known, several applications, in particular the analysis of phase 
transitions in the early universe and condensed matter physics require nonperturbative calculations. In this 
section we will derive an exact evolution equation for the dependence of the CGEA on the coarse graining 
scale [Q. In order to simplify the notation in this section we will denote by Aq the ultraviolet cutoff and 
by A the coarse graining scale. The CTP CGEA interpolates between the bare theory at A = Aq and the 
physical theory at the scale A. We will now consider a XcfA field theory in Minkowski spacetime. 



To derive such an exact evolution equation, we follow the approach of Wegner and Houghton 1 40 1 which 
is designed for Euclidean spacetime. Therefore we will invoke the Euclidean average effective action of 
Wetterich et al [ |37| , |67| -[7ll . (As already mentioned, the main difference between both actions is that the 
Euclidean action averages the field over a space-time volume, while our CTP CGEA averages the field over 
a spatial volume, and is therefore more adept to study non-equilibrium conditions.) 

The Euclidean CGEA is defined by 



g-SA(0) 



/ n P[0(g)] e-^-I-^l (VIII.l) 



Ao>(;>A 



By considering an infinitesimal variation A ^ A — (5 A it is possible to obtain an evolution equation for Sa, 
which reads 

,dSK ^ A f d^q f 5^Sa \ [' d^q SSaSSa f 5^Sa 
dA 26 A J (2^)4 i^<50,J0_ J J {27r)^ Scl)g 5q>.g yScl),6(ly^ J ^ 

(in the original equation of Wegner and Houghton there are additional terms coming from a rescaling of the 
modes after the coarse graining). The prime in the momenta integrals means that integration is restricted 
to the shell A> q> A - 5A 

A typical approximation taken to solve this evolution equation is to assume that is of the form 



Sa\4>] = / d^x 



(vin.3) 



Eq. ( VIII. 2 ) reduces in this case to an evolution equation for the effective potential Vi^{<j)). We will discuss 



this kind of approximation for the CTP CGEA in the next section. 
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Now we start the CTP calculation by writing the CGEA for a scale A — SA, namely 



ojSA-4A(0+,0-) — 



Ao>l«l>A-(5A 



(VIII.4) 



The modes to be integrated can be split into two parts: one within the shell A > \q\ > A — 5A and another 
containing modes with Aq > |g| > A. Expanding the action in powers of the modes within the shell, one 
obtains 



p'iS'A-5A('/'+,0-) 



giSA(</>+,0-) y. 



A>|g|>A-5A 



(VIII.5) 



where 



^1 

52 



= / dt 



(f'a{q,t) 



OS A 



- I dt dt' 



(27r)3 d(f)a{-q,t)(j)b[q,t') 



(VIII.6) 



In taking the functional derivatives of S'a (which contains modes whose wave vectors satisfy |(?| < A) the 
modes within the shell are set to zero. We use the notation 



4>a{q,t) 



Mq,t) 
<i>-{q.t) 



; 9ab — 



1 

-1 



(VIII.7) 



The S'a term is cubic in the modes within the shell and, as in the Euclidean case, it does not contribute 
in the limit JA ^ (basically, this is because one is doing a one loop calculation for the shell modes). The 
functional integrals over the she ll modes have the CTP boundary conditions. A comment about the last 
exponential factor in Eq.(VIII.5) is in order. Usually one discards it because it is a surface term, but in the 
CTP formalism it must be kept since the boundary conditions are that 0+ ((f, T) ~ 0_ ((f, T) with T — s- 00 
for the modes q within the shell. 

In order to evaluate the functional integrals we split the field as (f>a — (f'a + and impose the boundary 
conditions on the "classical" fields (j)±, i.e. they vanish in the past —T (negative and positive frequencies 
respectively) and match in the Cauchy surface at time T. The fluctuations vanish both in the past and 
in the future. The classical fields are solutions to 



(- 



dt^ 



q^)gab4'b{q,t) + / dt' 



where we have split the CGEA as SA{4>i 



dipai~q,t)dtpb{q,t') 

Skini(t>±) + Sint{4'±) with 



-M{q.t')^Q 



Ski' 



d'^X 



d^x 



(VIII.; 



(VIII.9) 



As before, in the functional derivati ves the modes within the shell are set to zero. 

Let ha be solutions to Eq.(VIII.S), vanishing in the past and satisfying an arbitrary normalization in the 
future, and let <p{q) be the common value of the fields taken in the future. We can then write 



4>a{q,t) = (l){q) 



haiq,t) 
ha{q,T) 



(VIII. 10) 



We first integrate over the common value (t){q) and then proceed with the functional integration over the 
fluctuations tpa (both are Gaussian integrals with "source" terms). One finally gets 
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A 



dSA 
dA 



iA 
2SA 



(2^ 



In 



h+{q,T) h_{q,T) 
h+{q,T) h.{q,T) 



A 
26A 



d^q h+{q,T) h^{q,T) 



{2^Y\h+{q,T) h^{q,T) 



iA ^ , I, , , 
— lndet'(A,b) 



A 



dt dt' 



d^q OSa 



dt 



(27r)3 dipa{q,t) 
The 2x2 matrix Aat has the foUowing elements 



ha{q,t) OSa 
' ha{q,T) dipa{-q,t) 



dSA 



(VIII.ll) 



A++{-q,t;q',t') = {-^ - q^ + ^e)5{t - t')S^iq + q') + 



dt 



A^^{-q,t-q\t') = {— + q^+te)6{t-t')6^q + q') + 



d^+{-q,t)d^+{q' ,t') 



A+_ (-g, t- q", t') = ', t'- -q, t) 



dip^{-q,t)d^^iq',t') 



dip+{-q,t)dip.{q'\t') 



(VIII.12) 



The primed determinant must be calculated as the product of the eigenvalues of Aab in a space of functions 
with wave vectors within the shell (A — (5A < |gj < A) and satisfying null conditions both in the past and in 
the future. Similar conditions are to be used to evaluate the inverse A~i^ . 



The equation (VIII.ll) is exact in the sense that no perturbative approximation has so far been used. It 



is similar to its Euclidean counterpart ( VIII. 2 ), but involves two fields and CTP boundary conditions. It 
contains all the information of the influence of the short wavelength modes on the long wavelength ones, and 
should be the starting point for a non-perturbative analysis of decoherence, dissipation, domain formation 
and out of equilibrium evolution. 



B. Derivative expansion 



The exact renormalization group equation is too complex to be solved without taking some approximation. 
The usual ones are expansions in the number of powers of the fields (see Ref. |Q for a detailed analysis) or 
in derivatives of them [f73|-[75[ . In the following we shall make use of th e derivative expansion approach. 

We will prove that, within this approach, the exact RG Eq. (VIII. 11) admits a solution of the form 



((/)+, (/;>_) = 5'a(0+) - Sa{4>-) 



(VIII. 13) 



Clearly this is not the most general form that can be imagined for the coarse grained action because contri- 
butions involving mixing of both fields are not taken into account. The main drawback of this approach is 
therefore that we miss the stochastic aspects of the theory, since there can be no noise terms. However, the 
proposed form for the CGEA will be enough for studying the renormalization group fiow of real time field 
theories. 



The great technical advantage of the form Eq. ( VIII. 13 ) is that the second functional derivative of the 
action has no crossed terms, leading to a diagonal matrix Aab whose determinant is easily computed as 

the product of two determinants, one for and one for A Following Ref. fz^ one can express both 

det' Aj^^ and det' A as the product over momenta of a constant (momenta independent) times the mode 

/i(g, T) evaluated at the final time T. Therefore the last term of the exact RG equation can be written as 



\ndet'{Aab) ^ ln[det'{A++)det'{A^^)] 



d\ 
(2^ 



■\n{h+(q,T)h_{q,T)) 



(VIII. 14) 



The first and the third terms can then be cast in the form of a single logarithm, and we arrive at 
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A 



dSA 
dA 




Hh^{q,T)h+{q,T) ~ h+{q,T)h^{q,T)) + 



h+{q,T) h_{q,T) 
h+{q,T) h^{q,T)^ 

(Pq OSa 



dt 



{27rr dipa{q,t) 



A-,\^q,t;q,t') 



ha{q,t) DSa 
ha{q,T) dipa{-q,t) 

dSA 



+ 



(VIII. 15) 



Note that the equations for the two modes and h_ (Eq. ( VIII.§| )) simphfy considerably, since the two 
equations are decoupled. What we still have to prove is that the proposed form for the action makes the 
r.h.s. of the exact RG equation split in the same form. 

Next we perform a derivative expansion of the interaction term. As our coarse graining explicitly breaks 
Lorentz invariance, we allow different coefficients for the temporal and spatial derivatives, namely 



4 



(VIII.16) 



We expand the fields around a time dependent background: <f>± = <j)± (t) + ip± (a;, t) and Fourier transform in 
space. We shall solve the Eq. (VIII. 8) for the modes to zeroth order in the inhomogeneities, i.e. we equate 
terms in the equations for h± that are independent of Vi's. Since the first functional derivative of the CGEA 
(S') is linear in the inhomogeneities ip±, we put S" = and keep the -independent contributions to S'/^f. 
After a little algebra and functional derivations, we get 



1 



dviq,t)d^{-q',t') 



5{t 



V" - ^<t> Z" - Yq^ - Z'i 
t')SHq-q') 



Z 



df^ 



4>z' 



(VIII. 17) 



where the primes denote derivation with respect to the field and the ellipsis denote terms linear in the fiuctu- 
ations. In this expression and hereafter we omit (unless explicitly stated) the ± subscripts in the background 
fields cj)±{t)^ in the potential VA(</)±(t)), and in the wave function factors ZA{<j)±{t)) and YA{(l)±{t))- Note 
that the effective mass of the modes depends on the time-dependent background (t>{t). The equations of 
motion for the modes ha become localized and take the form of harmonic oscillators with variable frequency 
and a damping term. The boundary conditions to be imposed are the aforementioned CTP ones. 

If one defines new modes as f{q, t) = {\ + ZAY^'^h{q, t), the damping terms cancel out and the new modes 
are harmonic oscillators with frequency 

1 + Ya 



^ ^ 1 + Za 1 + Za 

Using an adiabatic expansion for the modes, 

h±{q,t)^{l + Z,^-^/'- 



1 



7'2 



1 



Z'a 



4 (1 + Za)' 



21 + Za 



J'_^W±{q,t')dt' 



V^W±{q,t) 



(VIII. 18) 



(VIII.19) 



we can easily evaluate the logarithmic term in the r.h.s. of the exact RG E ( VIII. 15 ), which is the only term 
that survives in the approximation we are working. The RGE (VIII. 15) reduces to 



A / dt 



dVAih 



dA 



IdZAi^ 

2 dA ' 



dVAi^ 
dA 



i dZA{4>-: 

2 



dA 



A 
25A 



d^q 
(2^ 



[W+{q,t)-W-iq,t)]dt 



In the adiabatic expansion, the W's read 

=AA + BA4>''{t)+CAHt) 

where the coefficients are 



(VIII.20) 



(VIII.21) 
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Ba 



1 + Za 1 + Za 

1/2 
^A 



^A 



5A'? 



A" 
^A 



4(1 + Za)2 16v42 4Aa 



Ca- 



2(1 + Za) 4Aa 



Integrating by parts we get 

dVA 1 . dZ. 



dt I -A- 



dk 



2^ dk 



Ai2 




(VIII.22) 



(VIII.23) 



Therefore the dependence of the potential and the wave function renormalization on the infrared scale is 
given by 



A 



A 



dVA 
dA 

dZA 
dA 




(VIII.24) 



These equations are valid both for the 0+ field and for the 0_ field. 

The above equations describe the flow of the coarse grained action with the infrared scale in the derivative 
expansion of the exact CTP renormalization group equation. It is interesting to note that the higher 
derivative terms modify the differential equation for the effective potential. 

We have obtained two equations for the three independent unknown functions Va, Za and Ya- In order 
to find an additional relation between the spatial and temporal wave funtion renormalization functions Za 
and Ya, it is necessary to write the exact RGE up to quadratic order in the inhomogeneities. We will not 
present this long calculation here. For simplicity, we will assu me that Za and Ya sue small numbers, and 
therefore we will set them to zero on the r.h.s. of Eq.( VIII.24). This assumption is partially confirmed by 
numerical calculations . Note that in this approximation we recover the RG improved equation proposed 
in Ref. for the coarse grained effective potential 



A 



dVA 
dA 



A^ 

'47r2 



A2 



(VIII.25) 



There are other points which are worth noting. First, when we substitute Va, Za and Ya by their classical 
values Va = V , Za = Ya = 0, on the r.h.s. of both equations we obtain the one loop evolution equations 
p2t . Second, while in the one loop approximation it is possible to take the limit Aq ^ oo (the infinities can 
be absorbed into the bare mass and coupling constant), in this non-perturbative calculation it is not possible 
to renormalize the theory (as is the case for Hartree, Gaussian and 1/A^ approximations). For these reasons 
we keep Aq as a large (compared with the mass) but finite number. 

Once the functions Va and Za are known, one can write the effective dynamical equations for the coarse 
grained field. 



PART FOUR: Renormalization Group in Semiclassical Gravity 

IX. RENORMALIZATION GROUP AND STOCHASTIC SEMICLASSICAL GRAVITY 

In this Section we will describe a different relation between the RG and the CTP CGEA. We will consider 
the backrcaction of quantum matter fields on the spacetime geometry, assumed classical. To analyze this 
problem, we will use the formulation of quantum open systems and compute the CTP CGEA using an 
expansion in powers of the spacetime curvature and discuss its relation with the RG equations. 
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The usual approach to analyze backreaction in semiclassical gravity is based on the use of the Semiclassical 
Einstein Equations (SEEs) |£8[ 



1 



SttG 



- aHj,]) - PHj^J - ^ +<T^.>. (IX. 1) 



In the SEEs, which can be derived from the real part of the CTP CGEA , the effect of quantum matter 

fields is taken into account by including as a source the quantum mean value of the energy- momentum tensor. 
The terms proportional to 

ijW = [AR'"" - Agf^OR] + 0{R^), (IX.2) 

and 

hI^J = [4i?^"'^;„ - 2nRf"' - gf^'GR] + 0(i?2), (IX.3) 

come from terms quadratic in the curvature in the gravitational action, which are needed to renormalize the 
theory. 

These equations cannot provide a full description of the problem since they do not take into account 
the fluctuations of the energy momentum tensor around its mean value. The fluctuations can be incorporated 



by including an additional stochastic term [|7|, [79y80[ | on the right hand side of Eq. [X. 1 , This noise-term can 
be derived from the imaginary part of the CTP CGEA, in the same way we proceeded for the soft modes of 
the scalar fleld (see Section VII. B). When it is incorporated in the SEE, one obtains the "Einstein-Langevin 
Equations" (ELEs), which include both the dissipative and diffusive effects of the quantum matter on the 
geometry of spacetime [Q , in complete analogy with what happens in quantum Brownian motion typical of 
quantum open systems. For a review of the semiclassical stochastic gravity program based on the ELEs, see 

The ELEs have been derived for arbitrary small metric perturbations conformally coupled to a massless 
quantum scalar fleld in a spatially flat background ]80| ] , and, in a cosmological setting, for a massive field in 
a spatially fiat Friedmann- Robertson- Walker universe |^2|, and in a Bianchi type-I spacetime [Q. Further 
ellaboration of its physical meaning can be found in the papers by Verdaguer and Roura [ p4| . In ref. ]85| ] it 
is proven that the ELE may be used to compute certain quantum averages, even in conditions where there 
is no decoherence. 

Here we present the derivation of the ELEs based on the rcnormalization group equations of Lombardo 
and Mazzitelli p6[ |. Using a covariant expansion in powers of the curvature we will see that, for massless 
quantum fields, the ELEs are determined to leading order by the running couplings of the theory. 

This example will show an interesting connection between the "usual" rcnormalization group and the dif- 
ferent versions of the effective action. Indeed, in a naive "Wilsonian" approach, some quantum effects can be 
taken into account by replacing the coupling constants of the theory by their running counterparts. However, 
the (usual) running of coupling constants is defined in momentum space, and a careful implementation of the 
Wilsonian approach leads to a nonlocal effective action in configuration space. We will compute explicitly 
this non local effective action first in its Euclidean version, and show that it coincides with calculations 
based on resummations of the Schwinger DeWitt expansion. Using general relations between Euclidean 
propagators and CTP propagators, we will be able to obtain the CTP CGEA from the Euclidean effective 
action. From the CTP CGEA we will derive the ELEs, thus showing that they are a consequence of usual 
rcnormalization group equations of the theory. As an application, we will compute the leading quantum 
corrections to the Newtonian potential. 

Consider a quantum scalar field on a classical, Euclidean curved background. The classical action is given 

by 

S = Sgrav "I" Smatteri (IX. 4) 

where 

1 



Sgrav 



J d'^x^/g 
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{R ~ 2Ao) + aoi? + PoRt^.R'" 



(IX.5) 



and 
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(IX.6) 



Here ^ is the coupling to the curvature. Go, Aq, and the dimensionless constants ao and /3o are bare constants. 

The effective action for this theory is a compHcated, non local object. It is defined by integrating out the 
quantum scalar field, that is 

e"^"// j l?0e-'^lf^-"^l, (IX.7) 

where A'^ is a normalization constant. It is in general not possible to find a closed form for it. If we compute 
it using a covariant expansion in powers of the curvature, the different terms must be constructed with the 
Riemann tensor and its derivatives VV...7?.. Using integration by parts and the Gauss Bonnet identity, the 
effective action can be written only in terms of the Ricci tensor i?^^ and DR^j, |8^]. These arguments suggest 
that the effective action must have the general form 



S, 



eff 



R + aR^ + (3R,,^R'"' 



+ 



^ J d^x^[FoR + RFi{a)R + R^,F2{a)R^''' 



(IX.I 



where the ellipsis denote terms cubic in the curvature. For simplicity, in the above equat ion a nd in what 
follows we will omit the cosmological c onsta nt term. Note that the bare constants in Eq. ( [X.5| ) have been 



replaced by dressed couplings in Eq. ( |IX.8D . The expansion is adequate for weak gravitational fields, i.e. 
VV7^ < 7^2. 

Up to this order, all the information about the effect of the quantum field is encoded in the constant Fq and 
in the form factors Fi and F2. The form factors are, in general, non-local two point functions constructed 
with the d'Alambertian and the parameters ^ and m^. Fq, Fi, and F2 also depend on an energy scale /x, 
introduced by the regularization method. 

The dressed coupling constants depend on the energy scale /i according to the RG equations. 



minimal substraction these equations read 

dG G^m^ 



dfj, 

da 
^ dfjL 

dp 
^ dfi 



1 



327r2 
1 

960^ 



1 

90 



Using 



(IX.9) 



(IX. 10) 



(IX.U) 



The dependence of Fq, Fi and F2 on fi is such that the full equation is /i-independent. For example, from 
Eqs. (IIX.8D and ( |IX.9| ), we see that Fq = In (^^) - i) + const. 

When the scalar field is massless, this information is enough to fix completely the form factors. Indeed, as 
the Fi, i = 1,2 are dimensionless two point functions, by simple dimensional analysis we obtain Fi{0, jj? , ^) = 
Fi{^,£). Inserting this into Eq. ( IX.8| ), using Eqs. (IX. 10) and (IX.ll), and the fact that S^ff must be 
independent of /i, we obtain 



F,{U) 

F2ia) 
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const. 



const. 



(IX. 12) 



The final result for the effective action has a clear interpretation: it is just the classical action in which 
the coupling constants a and /3 have been replaced by nonlocal two point functions that take into account 
their running in configuration space. 

For a massive field, the situation is more complex because there is an additional dimensional parameter. 

2 

The form factors also depend on ^ and the /i-independence of the effective action is not enough to fix the 
form of them. They have already been computed in the literature |p8|,p9[ 
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where 



Xi(C,7) 



(IX. 13) 



^€(l-7') + 4(3-67^-7^) 



X2(^,7) = ^^7"*- 



(IX.14) 



These equations can be obtained through a covariant perturbation expansion |88[ , or by a resumma tion o f 
the Schwinger DeWitt expansion Of course these form factors coincide with our previous Eq. ( IX. 12 ) 
in the massless case. For the sake of completness we wih consider in what follows the massive case, although 

it is clear that only in the massless case the form factors are determined by the R G equa tions. 

In order to clarify the meaning of the two point functions appearing in Eqs. (IX. 12) and (IX. 13), it is 
useful to introduce the following integral representation 
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i(l-7')° 



In 



(1-7^) 



dz 
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(IX. 15) 



so the logarithm of the d'Alambertian is written in terms of the massive Euclidean propagator G^"* = 



4m^ 



□ ) ^. This representation will also be useful to construct the CTP version of the effective 



action. 

Up to here we considered the Euclidean effective action. What about the in-out and CTP CGEA? Of 
course one can compute them from first principles using the covariant expansion, and indeed there are some 
calculations in the literature for the CTP CGEA jo^. However, in order to emphazise the relation with the 
RG equations, we will construct the CTP CGEA from its Euclidean counterpart. 



Replacing the Euclidean propagator by the Feynman one in the integral representation Eq. ( IX. 15 ), one 
obtains the usual in-out effective action. As we already pointed out, the effective equations derived from this 
action are neither real nor causal because they are equations for in-out matrix elements and not for mean 
values. 

The CTP CGEA can be written as 



„.[s+,0+]-S„ 



.A9-,4>-]) 



(IX. 16) 



and the field equations are obtained from taking the variation of this action with respect to the metric, 
and then setting = g'^. 

In an alternative, and more concise notation, we can write this effective action as [ pl| 



49,0] 



(IX. 17) 



where we have introduced the CTP complex temporal path C, going from minus to plus infinity C+ and 
backwards C_, with a decreasing (infinitesimal) imaginary part. Time integration over the conto ur C i s 



defined by dt = /^^ dt — J^. dt. The field appearing in Eq. (IX. 17) is related to those in Eq. (IX. 16) 



by (/)(t, x) = (f>±{t, x) if t belongs to C±. The same applies to g^j^. 

This equation is useful because it has the structure of the usual in-out or the Euclidean effective action. 
Feynman rules are therefore the ordinary ones, replacing the Euclidean propagator by 



Gix,y) = 



Gf{x,v) = i{Q,in\T(l){x)(j){y)\Q,in) , 
GD[x,y) = -i{G,in\f(j){x)4>{y)\Q,in) 
G+{x,y) = -i{0,in\(j){x)(j){y)\0,in), 
G-{x,y) = i(O,m|0(y)0(a;)|O,zn), 



t, t' both on C+ 
t, t' both on C_ 
t on C_ , t' on C4 
t on C-I-, t' on C- 



(IX. 18) 



Introducing Riemann normal coordinates, we can write, up to lowest order in the curvature 



GF{x,y) 



(27r)4 p2 



= G*D{x,y), 



(IX. 19) 
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(IX.20) 



All of the preceeding formulation of the effective action is valid for any field theory. In our particular case, 
we must replace the Euclidean propagator in Eq. (IX. 15) by the propagator G{x,y) of Eq. (IX. IS) 
with a mass given by ^™ ^ + z. After integration in z we obtain 



In 



4m^ 



CTP 



f APpip{x-V) It, ( (l-7")(p"7»e)+4m^ 



4m 



i, i' both on C+ 
i' both on C_ 
i on C_ , i' on C_|_ 
t on C+, t' on 



(IX.21) 



With the expression for the CTP logarithm of the d'Alambertian we can calculate explicitly the CTP 
effective action. Using the previous notation with and the CTP effective action reads 



'8^ 



d^x / d''yA{x)A{y)Ni{x,y) 



1 

8^ 



d^x / d''yA{x)J:{y)Di{x,y) 



J d'x j dSA^,ix)A^^y)N2ix,y) - ^ J d^x J d^yA.^ix)^'^'' iy)D2{x,y), (IX.22) 



where A = 



-R7, 



2 , ^ — 2 ' "^A"' ~ '"'2 J '^p.y = '^"2 ■ "^^^ classical gravitational action Sg^^^ 
contains the dressed, /x-dependent coupling constants and we absorbed Fq into the gravitational constant G. 

The real and imaginary parts of Seff can be associated with the dissipation and noise, respectively. The 
dissipation Di and noise Ni kernels are given by 



Dz{x,y) = / d'^Xii^n) 



_dS_ 
(2^ 



cos[p(x — y)] In 



(1 - 7^)^^ + 
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N.i{x,y) 



cos[p(x — y)\{) —p 
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1-7^ 
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It is important to note that the imaginary part of this effective action must be positive definite. To make 
this point explicit, one can write the imaginary part in terms of the Weyl tensor G^^ap and the scalar 
curvature R by means of the following relation: C^^a/jC'"'"'^ = 2R^,^R^"' - 2/3i?2. It is not difficult to show 
that the scalar and tensor contributions to the imaginary part of the effective action are both positive. 

In order to derive the ELEs we proceed as in Section VII B. One can regard the imaginary part of the 
CTP CGEA as coming from two classical stochastic sources ri{x) and ■q^'^"'^{x), where the last tensor has 
the symmetries of the Weyl tensor. In fact, we can write the imaginary part as 



V7^{x) / V,^^^''^P{x)P[7^, rj'^'"'^] exp {i {A(a;)7?(x) + A,,.,^??''''"^}) 



exp < — d X d y 



'a{x)N{x - y)A{y) + A^,^p{x)N2{x - y)A'^^''^{y) 



(IX.25) 



where N{x,y) = Ni{x,y) + ^N2{x,y), and A^^^p — ^C'^^^^ — G^^^^p. The Gaussian functional probability 
distribution P[ry, jy^'^"^] is given by 



P[V, »r"^] = A exp (^-^ J d^x J d^yr^ix) [n{x, y)] ' rj{y) 

xexp|-i / d'^x I d^yri^^^pix) [N^ix.yT^^^^^ {y) 



(IX.26) 
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with A a normalization factor. 

Therefore, the CTP CGEA can be written as 



where 



exp{iS'e//} = J I??7l??7^^a/3P[77, 77^^0/3] exp{iAe//[A,A^i,„^,i;,i;^i,,?7,77^i,„^]}, (IX.27) 
Aeff^ReSeff+ [ d^x[A{x)7j{x) + A^,^0{x)t]^"'"'^]. (IX.28) 



The field equations 



if 



= 0, the Einstein-Langevin equations, are 
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dSD2{x,y)H^^J{y) 



2r] 



a/3 



(IX.29) 



where a and (3 differ from a and (3 by ^-dependent finite constants. Eq. ([X.29) is the main result of this 
section. The r.h.s. consists of the mean value of the energy-momentum tensor for the scalar field plus a 
stochastic correction characterized by the two point correlation functions 



< r]p.:^a(3{x)r]p^XTiy) > 



N{x,y) 

Tfj.ua!3pcrXTN2{x, y) , 



(IX.30) 



where the tensor T, 



^ is a linear combination of four-metric products in such a way that the r.h.s of 



Eq. ( IX.30 ) keeps the Weyl's symmetries. The scalar-noise kernel is given by 



N{x,y) - - 



^7 



(1-7^) 
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In the massless case N is proportional to — 1/6)^, and vanishes for conformal coupling. Therefore this 
term is present when the quantum fields are massive and/or when the coupling is not conformal. This 
is to be expected, since the imaginary part of the CTP CGEA signifies particle creation. For massless, 
conformally coupled quantum fields, particle creation takes place only when the spacetime is not conformally 
flat. Therefore in this case the only contribution to the imaginary part of the CTP CGEA is proportional to 
the square of the Weyl tensor. When the fields are massive and/or non-conformally coupled, particle creation 
takes place even when the Weyl tensor vanishes. This is why an additional contribution proportional to i?^ 
appears in the imaginary part of the effective action. 

From Eq. |X^) we can define the effective energy-momentum tensor 



rpstoch 



=<T^, >+g'^'^□7/-7^• 



,;Mf 
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(IX.32) 



where < T^jy > is the quantum expectation value of the energy- momentum tensor of the quantum field 
and T^*°'^^ is the contribution of the stochastic force, which in turn has contributions from the scalar and 
tensor noises. In the massless-conformal case the scalar-noise kernel vanishes, and {T^^^y^o<^h _ because 
the noise-source r]^'^°'^ has vanishing trace. This means that there is no stochastic correction to the trace 
anomaly |^ . 

To summarize, we have obtained the ELEs using a covariant expansion in powers of the curvature. Our 
results are valid for quantum scalar fields with arbitrary mass and coupling to the curvature ^. In the massless 
case, still for arbitrary ^, we have shown that it is possible to obtain the noise and dissipation kernels using 
only dimensional analysis, the running of the coupling constants and the relation between the Euclidean 
and CTP CGEA. From this point of view, we can therefore conclude that the RG equations already contain 
information about the dissipation and noise kernels. 

The results of this section can be easily generalized to any field theory with massless quantum fields in 
gravitational or Yang Mills backgrounds. 
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X. RENORMALIZATION GROUP AND QUANTUM CORRECTIONS TO THE NEWTONIAN 

POTENTIAL 



Two major areas of application for the stochastic and semiclassical gravity programs are the study of the 
final states of black holes and the physics of the early Universe. These are enormously complicated problems. 
We will describe here one of the simplest applications of semiclassical gravity, the calculation of quantum 
corrections to the Newtonian potential. Although a very simple example, it will shed light on important 
conceptual issues, in particular on the relation between the quantum corrections and the renormalization 
group. 

Duff computed the corrections to the Newtonian potential produced by the vacuum polarization of 
gravitons. His result is, schematically, 



Veffir) 



GM 



(1 



G 



(X.1) 



This result has been rederived more recently by a number of authors ]93|-|95||. It was also shown that 
corrections of this form are the leading quantum corrections for massless quantum fields. The constant a 
depends on the number and spin of massless fields. The first derivation based in the SEE has been presented 
in Ref. ]9^ . Recently there is a renewed interest in this type of corrections since it is relevant to relating 
two different developments in quantum gravity: Maldacena's Ads/CFT correspondence and the Randall- 
Sundrum alternative to compactification, see for example [ p7| . 

There is an alternative, intuitive "Wilsonian" way of taking into account, at least partially, the quantum 
effects: just modify the classical potential by replacing the Newton constant by its running counterpart 



V{r) 



(X.2) 



where G{n) is th e solu tion to the renormalization group equations in the theory considered. As can be seen 
from Eqs. ( |IX.E ) and this argument does not reproduce the leading quantum correction for massless 

fields. 



Wc will now show that the result ( X.1 ) can be derived from the ELEs, and are therefore a consequence of 
the RG equations for the parameters a and (3. For simplicity, and to make contact with previous works, we 
will solve the ELEs without including the noise source. Solutions to the ELEs including this source can be 
found in g. 

In the static, weak field approximation g^u — -q^j^i, + ft,^^, we have 
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^■^ fJ,V 
^(2) 



-nh 
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(-29^9^ + 277^^D)n/i 



where h = rj'^'^hf^^ and we assumed the Lorentz gauge condition (/i^^ — hhri^i^y = 



Including a point particle source with T^^ = J^5"(5'^(x), the linearized field equations become, 
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(X.4) 
(X.5) 
(X.6) 
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with hfj_i, = hfj_i, — i^h-q^y. We define the quantum corrected Newtonian potential as V{r) 
The trace of the field equations is 
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where, to first order in 
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(X.8) 
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We shall solve Eq. (X.8) perturbatively h — /i^^^ + h^^^ . The classical contribution satisfies 
(V2 - r7-2v2v2)/i(") = -16^GAf(53(x) ^ 327rG(3a + /3). 

The time independent and spherically symmetric solution is 

4GM 
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The equation for the first quantum correction is 
where 
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To find a solution to this equation, we will consider the limit or — )■ oo (we are interested in long distance 
quantum corrections). In this limit h^^^ = 4GM + 4:TTa~^S^{x)) and the nonlocal operator in the rhs of 
Eq.( X.l^ can be easily evaluated. After a long calculation we obtain 
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The dots denote corrections at the origin, which are proportional to (5"^(x) and its derivatives. We have not 
included them because our quantum corrections are not accurate near the origin. Indeed, we have derived 
the modified Einstein equations under the assumptions WTZ ^ TZ"^ and m^TZ ^ VVi?. Both conditions 
are satisfied for the potential if GM <C r- ^ to~^, so the origin is excluded. 

A similar analysis can be carried out for hoQ. We omit the details |Q. The final answer for the Newtonian 
potential is 



V{r) - ~-hoo = 
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(X.15) 



From Eq. (X.15) we see that there are two different terms in the quantum correction. The term containing 
the logarithm is qualitatively what we expected from 'Wilsonian' arguments (sec Eq. X.2). Howe ver, the 
coefficient is not exactly the same as the one derived from the renormalization group equation ( IX.9| ), unless 
^ = 0. Besides the running of G, we have obtained the leading quantum correction, which is a consequence 
of the running of a and (3 in configuration space. 



XI. RENORMALIZATION GROUP THEORY FOR NONEQUILIBRIUM SYSTEMS 

A. Summary Remarks 

We first give a schematic summary of this report and then discuss some general issues concerning RG 
theory for NEq systems. 

The CGEA contains all the information about the influence of the environment on the system. The CTP 
version of the CGEA is suitable for the analysis of the dynamical evolution of the system under the influence 
of the environment, and takes into account effects of renormalization, dissipation and noise. 

In this review we have described several applications of the in-out and in-in effective action in semiclassical 
gravity and cosmology, paying particular attention to the relation with the RG equations. The Euclidean 
averaged effective action is reviewed by Wetterich and his collaborators in an article in this conference. 

We have described a perturbative calculation of the CTP CGEA in RW spacetimes. We have shown 
that the CTP CGEA is useful to derive the master and Langevin equations in quantum field theory. The 
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imaginary part of the CTP CGEA produces diffusive terms in the master equation. These terms can induce 
the reduced density matrix to become diagonal as it evolves, and are therefore crucial for decoherence and 
quantum to classical transition of the system. We applied this formalism to investigate the decoherence of 
mean fields due to quantum fluctuations relevant for theories of structure formation in inflationary models. 

We adapted the Wegner and Houghton Euclidean approach to the in-in CGEA in order to obtain an 
exact renormalization group equation for the dependence of Sa on the coarse graining scale. The exact 
equation is extremely complicated. We have solved it using a derivative expansion. In this approximation, 
the CTP CGEA "decouples" [i.e. it is of the form Sa{4'+, 4>-) = 5'a(^+) — S\{4>-) ] and contains neither 
dissipative nor noise terms. Previous results on the RG improved effective potential are recovered under this 
approximation. 

Efforts to find solutions beyond the adiabatic approximation are under investigation. We believe that this 
equation (or a simplified version of it) will play an essential role in the development of a renormalization group 
theory for nonequilibrium systems. We expect that, as soon as we decrease the scale from Ag, dissipative 
and noise terms will grow: the CGEA will develop an imaginary part (related to noise) and a real part 
containing interactions between the (j>± fields (dissipation). This can be easily checked both in the one loop 
approximation and from the exact RGE. Indeed, we have seen that the one loop CTP CGEA is i n genera l 
non-real. On the other hand, the real and imaginary parts of the CGEA are not decoupled in Eq. ( VIII. 11 ), 
and a non- vanishing real part at A = Aq will induce an imaginary part at lower scales. 

One should be able to find a non-perturbative, A-dependent fluctuation-dissipation relation and RG equa- 
tions for the coupling constants of the theory that include the noise effects. 

Finally, in the context of semiclassical gravity we described an interesting relation between the usual 
RG equations (running coupling constants) and the CTP CGEA. We have seen that one can take into 
account backreaction effects of quantum fields on the spacetime metric using a "Wilsonian" effective action 
in which the parameters of the theory are replaced by their running counterparts (this form of the effective 
action suggests itself by demanding the theory to be independent of the scale /i introduced by dimensional 
regularization). As the running is local in momentum space, it becomes nonlocal in configuration space, 
and thus one obtains an Euclidean, nonlocal effective action for the spacetime metric previously obtained 
by using different approaches. This Euclidean action can be transformed into the CTP CGEA taking into 
account formal analogies between their definitions. The final result, which contains dissipation and noise 
effects, can be viewed as a consequence of the running couplings of the theory. 



B. Towards a Nonequilibrium Renormalization Group Theory 

There are three important aspects in the construction of a renormalization group theory for nonequilib- 
rium systems: a) The idea behind the introduction of the renormalization group (RG) description of critical 
phenomena characterized by the running of the interaction parameters of the theory; b) The scale charac- 
teristic of the dynamical interaction absent in static critical phenomena; c) How noise and dissipation in the 
stochastic equations describing the effective dynamics are reflected in the renormalization group equations 
governing the parameters of the theory. Point a) is discussed in all theories of static critical phenomena. 
We approached this problem stressing t he co arse graining and backreaction aspects. Point b) is discussed in 
dynamical critical phenomena (see e.g., |lOCf| ). We note, however, in the literature the noise in the dynamical 
Ginzburg-Landau equation is usually introduced by hand, and what is at issue is how the system behaves 
towards the critical point ~ neither is the origin of noise accounted for nor is the dissipative aspect of the 
system dynamics being incorporated in the flow. The first part of Point c), i.e., how dissipation and noise 
appear in an open system is discussed in this report. (For a discussion of how noise is identified in interacting 
quantum fields in both a prescribed open system and an effectively open system, see, |8 
remaining issues are how they affect the RG fiow and how they manifest in the approach to critical points. 
We will make some general comments here, as research on nonequilibrium (NEq) RG is still in its infant 
stage, its application to self-organized criticality, driven-diffusive systems and turbulence notwithstanding. 

It is perhaps helpful to refiect upon the basic ideas and procedures behind the introduction of the RG - 
i.e., coarse graining and scaling. When a system has a certain degree of regularity (homogeneity, periodicity) 
depicted by some symmetry, one can choose to represent it with an equivalent coarse-grained description. 
Examples are decimation in real time RG applied to a lattice, where, e.g., in the Ising model every other spin 
is eliminated and the effective bond strength between remaining spins is doubled. This Kadanoff-Migdal 
procedure produces a new Hamiltonian, which, upon iteration, can transform an original system into a simpler 
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one. In the momentum space description, this procedure transcribes (runs) the ultraviolet (short range) 
behavior of the system to its infrared (long-ranged) domain, and hence is useful for the description of critical 
phenomena, as the behavior of the system near the critical point is dominated by the appearance of long range 
order. Whether this procedure, which is one of many possible coarse-graining schemes, can faithfully depict 
or capture the essential physics depends on the compatibility of the procedure with the system and on the 
properties of the system near the critical point. If the coarse- graining respects the symmetry of the system, 
(e.g., for an anisotropic medium, use a different decimation grading in different directions commensurate 
with the symmetry of the medium) and if the system possesses some scaling properties near the critical 
point, then the transformed problem via the RG would preserve the same critical behavior as the original 
problem. Otherwise it fails. The applicability of the RG idea thus depends crucially on the choice of a 
coarse-graining procedure and the use of scaling concepts. Leaving scaling aside for now, which has more to 
do with the properties of the particular systems of interest than with the procedure for accessing relevant 
information about the system, our approach to nonequilibrium processes and NEqRG theory starts from 
examining closely the coarse-graining procedure. 

1. RG procedures in the light of open system concepts 

As we stated clearly at the beginning of this report, the necessary steps to capture the essense of a physical 
system with a simplified depiction lies in : 1) distinguish the system from the environment; 2) coarse grain 
the environment and 3) measure how the coarse-grained environment influences the system in providing an 
effective kinematics or dynamics of the reduced system. Let us examine the RG procedure in the light of 
the open systems scheme of nonequilibrium statistical mechanics. The first step consists of a) separating the 
order parameter field ^{x,t) into two parts, $ — <I>s + ^e, where $5 and are respectively our system 
and environment. We assume that $s contains the lower k wave modes and ^e the higher k modes. In 
critical phenomena, -l k |< A/s, ^e '■ A/s <| k |< A. Here A is the ultraviolet cutoff and s > 1 is 
the coarse-graining parameter which gives the fraction of total k modes counted in the environment. Then 
b) one 'integrates out' the short wavelength sector, which amounts to finding out the effect or backreaction 
of the coarse-grained environment (the short wavelength modes) on the system. This is most succinctly 
and forcefully executed by the use of coarse-grained effective action as illustrated in the examples of 
this report. The last step in the RG procedure is c) to introduce a rescaling of the k-space and a rescaled 
field, e.g. k' = sk, (j>'g{k',t) — s ~(j)s{k,t) (d being the spatial dimension). In identifying (or rather, 
demanding that they are equivalent, up to a certain order in the perturbation expansion) the new action in 
terms of the rescaled variables with the old action, one can obtain a set of differential RG, or the Wegner- 
Houghton equation. At the critical point itself, this effective theory containing only the low frequency modes 
is insensitive to further coarse graining, and therefore is described by a fixed point of the RG flow. 

For dynamical systems, in addition to the scaling transformation, an artificial device of the RG theory 
depicting the running of the coupling constants of the system towards the infrared region, we find a real 
time-dependence depicting the dynamics of the system. ^ 

Now that we have identified the steps of the RG procedure in parallel to the treatment of open systems we 
can ask a few questions to illuminate the relevant issues. It is easy to see where dissipation and noise occur 
in the open system. But the appearance or nonappearance of noise and dissipation in the RG equations for 
nonequilibrium dynamics is not a simple issue. That is why we deem it necessary to first provide a clear and 
thorough discussion of the conceptual and technical basis for addressing the origin and nature of noise and 
dissipation in nonequilibrium field theory, as this report attempts to do, before we move on to investigate the 
related issues in RG theory for nonequilibrium processes. When the RG procedure is viewed in the light of an 
open system theory — in effect the modes of interest in the system after the Kadanoff transformation form 
an open system, and the techniques of nonequilibrium field theory are the natural language to address these 
problems — the long wavelength / slow modes dynamics will display, in general, both noise and dissipation. 



In treating dynamical systems, we may also coarse grain the system according to some time scale, obtaining the 
effective theory for slow modes after coarse graining the fast modes. Sometimes what is slow compared to what is 
fast may not be taken at face value, but has to be determined relevant to the dynamics and symmetry of the system. 
(Viewing eternal inflation as 'static' and 'slow- roll' as dynamic is an example we showed.) 
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Indeed, the possibility that a renormahzation group flow may acquire stochastic features for nonequilibrium 
processes is even clearer if we think of th e RG as encoding the process of eliminating irrelevant degrees of 
freedom from our description of a system |102|. These elimination processes lead as a rule to dissipation and 
noise, as was made manifest in the influence action and the CTP-CGEA approaches shown. 

But why, one may ask, is it that we don't usually talk about dissipation in the system or noise in the 
environment in a RG theory. A simple answer is that RG running is different from real-time dynamics, and 
dissipation in the dynamics of the effective system does not show up as dissipation runs. Another simple 
answer is that the bulk of R G res earch has been focused on equilibrium, stationary properties rather than 
the nonequilibrium dynamics [ |l03[ . f\ A major task in NEqRG is to make explicit this point, i.e., where does 
dissipation in the open system (Langevin) dynamics show up in the RG equations? We can offer only some 
speculative observations: 



2. Stochastic RG Equations 

A) For apparently closed (yet effectively open) systems described by the Boltzmann dynamics, this issue 
is more involved. We learned that when the hierarchy of correlation functions are simply truncated (with no 
slaving) |101| ] , the equation of motion for the finite set of low-order correlation functions is unitary, such as is 
the case in Vlasov dynamics, which is obviously non-equilibrium, yet no dissipation or noise appear as such. A 
similar case is the n-loop effective action in ordinary quantum field theory. If we view the separation between 
the classical background field and the quantum fluctuation fleld as between a system and an environment, 
then the ordinary n-loop effective action is only one special case of the coarse-grained effective action, with 
quantum fluctuations being integrated or coarse-grained away. This results in radiative corrections to the 
bare mass and charge as we learned from renormahzation theory, which, in the nonequilibrium language, 
amounts to backreaction of the environment on the system. We would expect to see dissipation in the effective 
dynamics of the background flcld if we follow the general arguments of NEq fleld theory, but of course one 
never talks about dissipation or noise in the equations of motion derived from an n-loop effective action. So 
what is missing? The reason is similar for both cases: the mean fleld or background field is the lowest order 
of the hierarchy of correlation functions in the Schwinger-Dyson equations. Any loop expansion entails a 
truncation of the hierarchy and, just like Vlasov dynamics in classical mechanics, the equation of motion 
from the effective action is unitary. So again it is not just the dynamics which prompts the appearance of 
dissipation - the causal factorization condition which wc call slaving (as in the molecular chaos assumption 
of Boltzmann) is responsible for it. 

For example the circumstances whereby the ordinary RG equations are derived could be similar to the 
Vlasov or the loop approximation examples mentioned above, in that the correlation functions are trun- 
cated from the hierarchy with a factorization condition (rather than from slaving, where the correlation 
noise arises). This results in the coupling constants being modified only partially, in a way similar to the 
renormalized mass or charge in one-loop (quantum) fleld theory or the (classical) Vlasov dynamics following 
the averaged potential which replaces the particle interactions. In other words, the backreaction of the 
coarse-grained modes is not taken into account fully. We see this, e.g. for the tj)"^ theory: the rescaled theory 
differs from the original theory in the 0^ and higher order terms. By identifying the two theories in such a 
way prescribed by the ordinary RG theory procedure, we assumed that the error of ignoring the higher order 
corrections are inmaterial in our range of consideration. As such it is a good description of the theory only 
for systems which scale near the critical point, because there the higher order correlation functions can be 
ignored or incorporated in a way simply related to that of the lower order ones. Away from the critical point 
this is not true and this procedure would give bad results. That is the criterion whereby one can assert the 
equivalence of the RG transformed theory and the original one. 

Thinking about the relation between the general nonequilibrium dynamics of an open system and that 
depicted by the RG transformation in this light, we can see that the former (say, starting from a nonlocal 
Langevin equation governing the system arising from the backreaction of another subsystem or environment 



® As for whether it has to do with equilibrium versus nonequilibrium conditions, we know that even under equilibrium 
conditions there is always noise and dissipation, as can be seen from linear reponse theory. They are related by the 
fluctuation-dissipation theorem. Indeed it is the precise balance of these processes which sustains the equilibrium. 
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via a Zwangzig-Mori projection operator) keeps track of much more information about the environment than 
that in the RG treatment of the critical regime (as slaving being more than truncation). Although this is not 
sufficient for a proper description of the effective dynamics as dissipation and fluctuations carry important 
information about the system and the environment, it is sufficient for the purpose of depicting the (static) 
critical phenomena of the system in the critical region, as scaling behavior simplifies the relation between 
the correlation functions. That is why in actual RG calculations one can assume that the cutoff is in effect 
infinitely far above the scales of interest. This advantage near the critical point will still come to the aid 
of dynamical critical phenomena, but when the full dynamics goes beyond an approach to the equilibrium 
(which is manifest in say, critical slowing down) or when one is interested in the behavior of the system 
farther away from the critical point, then more of the dissipative and noise attributes which are always 
present in general circumstances will have to be included in a nonequilibrium RG theory. The art in the 
design of such a theory for NEq processes is to find out what additional information (noise and dissipation) 
need be kept, or what degree of coarse-graining is sufficient for the RG procedure which can still capture the 
essence of the critical behavior for dynamical systems. 

B) One place where one can indeed see the 'disappearance' of the effects of noise from the RG equations 
is when an adiabatic or quasilocal approximation is introduced on the mode function dynamics. In 
Dalvit and Mazzitelli show from the Wegener - Houghton equation derived from the CTP effective action 
that, for static long wavelength fields, the usual RG is recovered. In view of this result it is clear that an 
understanding of the running of noise and dissipation requires consideration of dynamical long wavelength 
fields. However, if these fields are slowly varying, they may be considered as actually linear or quadratic in 
the space - time coordinates, in which case the fu nctio nal integration of the high momentum modes may be 
carried out exactly, at least to one loop order |p6,104|. 

Alternatively, we may consider that the low momentum modes are following a prescribed, though not 
necessarily slow, evolution, and compute the back reaction of the high momentum modes keeping track 
of non adiabatic effects. Such an analysis is carried out in [f05|(See also Q), where the fast modes are 
represented by a massive heavy field interacting with slow modes, described as a light field. As expected, 
noise and dissipation display significant running (they increase exponentially) , as the gap between the char- 
acteristic frequencies of the heavy and light fields narrows. This result suggests that a nonequilibrium RG 
should include, besides the usual interactions, the running of suitable parameters describing both noise and 
dissipation. 

Therefore the challenge is to construct a RG flow where the dissipative and noisy characteristics of the 
dynamics is ingrained therein. There is a precedent for this problem in studies of RG flow in the Navier 
- Stokes equation [ 106 |, where one is interested, for example, in the running of the viscosity. In this kind 
of research, the starting point is generally a Langevin type equation , wh ose solutio ns a re given through 
functional representations following the methods of Zinn- Justin, e.g., 100 (see also [107|). We would like 
to point out that even in this context, closed time - path techniq ues o ffer an alternative way to obtain a 
functional representation of the solutions of the Langevin equation |108|, which is in several aspects simpler 
than the more familiar, forward time path, one. 

C) Another lead can be found in the stochastic correlation function G introduced in a quantum field 
theoretical derivation of the Boltzmann-Langevin equation |101]. Its expectation value reproduces the usual 
propagators (Green functions), while its fluctuations account for the quantum fluctuations in the binary 
product of (operator) fields. The dynamical equation for G takes the form of an explicitly stochastic Dyson 
equation. In the kinetic limit, the fluctuations in G become the classical fluctuations in the one particle 
distribution function, and the dynamical equation for G's Wigner transform becomes the Boltzmann - 
Langevin equation. (Each of these results has an interest of its own. A priori, there is no simple reason 
why the fluctuations derived from quantum field theory should have a physical meaning corresponding to 
a phenomenological entropy flux and Einstein's relation.) The studies of the fluctuating character of these 
field theoretic Green functions also suggest new avenues in the development of RG theory. For example, we 
are used to fixing the ambiguities of renormalization theory by demanding certain Green functions to take 
on given values under certain conditions (conditions which should resemble the physical situation of interest 
as much as possible, as stressed by O'Connor and Stephens |10£]). If the Green functions themselves are 
to be regarded as fiuctuating, then the same ought to hold for the renormalized coupling constants defined 
from them, and for the renormalization group (RG) equations describing their scale dependence. The notion 
that Green functions (and indeed, higher correlations as well) may or even ought to be seen as possessing 
fluctuating characters (when placed in the larger context of the whole hierarchy) with clearly discernable 
physical meanings is likely to have an impact on the way we perceive the statistical properties of field theory. 
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To end, we note that while the application of rcnormalization group methods to stochastic equations is 



presented in well-known monographs |10C], our proposal here goes beyond these results in at least two ways. 
First, in our approach the noise is not put in by hand or brought in from outside (e.g., the environment of an 
open system), as in the usual Langevin equation approach, but it follows from the (quantum) dynamics of 
the system itself. Actually, the possibility of learning about the system from the noise properties (whether 
it is white or coloured, additive or multiplicative, etc.) - unraveling the noise, or treating noise creatively- 
is a subtext in our program. Second, our result suggests that stochasticity may, or should, not only appear 
at the level of equations of motion, but also the level of the RG equations, as they describe the running of 
'constants' which are themselves fluctuating. This we feel is the flrst task in the construction of a RG theory 
for NEq processes. 
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